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Abstract  (coat) 

This  Semiannual  Report  also  presents  statistics  from  operation  of  the  Intelligent  Monitor¬ 
ing  System  (IMS).  The  IMS  has  been  operated  in  an  experimental  mode,  and  the  perfor¬ 
mance  has  been  very  satisfactory.  Since  October  1991,  a  new  version  of  the  IMS  that 
accepts  data  from  an  arbitrary  number  of  arrays  and  single  3-component  stations  has  been 
operated. 

The  NORSAR  Detection  Processing  system  has  been  operated  throughout  the  period  with 
an  average  uptime  of  97.3%  as  compared  to  97.8%  for  the  previous  reporting  period.  A 
total  of  2290  seismic  events  have  been  repealed  in  the  NORSAR  monthly  seismic  bulletin. 
The  performance  of  the  continuous  alarm  system  and  the  automatic  bulletin  transfer  by 
telex  to  AFTAC  has  been  satisfactory.  The  system  for  direct  retrieval  of  NORSAR  wave¬ 
form  data  through  an  X.25  connection  has  been  used  successfully  for  acquiring  such  data 
by  AFTAC.  Processing  of  requests  for  full  NORSAR  and  regional  array  data  on  magnetic 
tapes  has  progressed  according  to  established  schedules.  There  have  been  no  modifica¬ 
tions  made  to  the  NORSAR  data  acquisition  system. 

On-line  detection  processing  and  data  recording  at  the  NORSAR  Data  Processing  Center 
(NDPC)  of  NORESS,  ARCESS  and  GERESS  data  have  been  conducted  throughout  the 
period.  FINESA  data  were  processed  until  16  May,  when  field  work  for  the  array  refur¬ 
bishment  started.  Data  from  two  experimental  small-aperture  arrays  at  sites  in  Spitsbergen 
and  Apatity,  Kola  Peninsula,  have  been  recorded  and  processed  in  an  experimental  mode. 
Monthly  processing  statistics  for  the  arrays  as  well  as  results  of  the  IMS  analysis  for  die 
reporting  period  are  given. 

Maintenance  activities  in  the  period  comprise  preventive/corrective  maintenance  in  con¬ 
nection  with  all  the  NORSAR  subarrays,  NORESS  and  ARCESS.  Other  activities  have 
involved  testing  of  the  NORSAR  communications  systems  and  work  in  connection  with 
the  experimental  small-aperture  arrays  in  Spitsbergen  and  Russia. 

Starting  1  October  1991,  an  effort  began  to  carry  out  a  complete  technical  refurbishment 
of  the  NORSAR  array.  This  project  is  funded  jointly  by  AFTAC,  ARPA  and  NFR.  During 
the  reporting  period,  we  have  continued  evaluation  and  laboratory  testing  of  technical 
options  for  field  instrumentation,  in  particular  state-of-the-art  A/D  converters,  data  acqui¬ 
sition  and  synchronization  devices.  During  the  reporting  period,  we  have  also  received  a 
complete  subarray  evaluation  unit  from  Science  Horizons,  and  this  unit  has  been  subjected 
to  a  long-tom  operational  test  in  subarray  06C.  In  addition,  we  have  been  converting  the 
NORSAR  Detection/Event  Processing  software  to  be  compatible  with  our  UNIX-based 
workstations.  Our  detailed  status  regarding  this  work  is  reported  on  separately,  most 
recently  in  a  status  report  dated  16  August  1993. 

Summaries  of  six  scientific  contributions  are  presented  in  Chapter  7  of  this  report. 

Section  7.1  describes  an  experiment  in  continuous  threshold  monitoring  of  the  Lop  Nor, 
China,  nuclear  test  site.  The  monitoring  period  comprised  the  five  days  up  to  and  includ¬ 
ing  the  day  of  the  latest  nuclear  explosion  at  Lop  Nor  (5  October  1993).  NORESS, 
ARCESS  and  GERESS  regional  array  data  were  used.  We  have  been  compiling  daily  sta- 
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tistics  of  til  peaks  on  the  threshold  diagram  exceeding  mb  3.S  and  associating  these  peaks 
to  regional  or  teleseismic  events  whenever  possible.  Based  on  this  experiment,  it  appears 
that  this  threshold  (3.5),  which  is  about  0.75  0%  units  higher  than  for  Novaya  Zemlya,  is 
representative  far  the  combined  threshold  monitoring  capabilities  (at  the  99  per  cent  level) 
of  these  three  regional  arrays. 

Section  7.2  presents  some  observations  for  the  mb  =  5.9  Lop  Nor  explosion  on  5  October 
1993.  This  explosion  was  naturally  observed  on  the  threshold  plot  of  Section  7. 1  as  a  very 
sharp  peak.  The  contribution  summarizes  the  automatic  processing  results  for  all  six 
arrays  processed  at  Kjeller.  The  NORESS  array  has  the  best  signal-to-noise  ratio  (1376.2) 
for  this  event,  and  by  extrapolation  this  array  would  be  expected  to  have  a  detectable  sig¬ 
nal  for  an  event  from  Lop  Nor  about  2.5  magnitude  units  lower.  ARCESS  and  NORSAR 
also  show  outstanding  SNR.  The  velocity/azimuth  estimates  are  within  the  expected 
uncertainty  for  all  arrays.  It  is  also  noted  that  a  large  earthquake  and  aftershock  sequence 
occurred  in  the  Southern  Sinkiang  province  only  three  days  before  the  explosion,  and  this 
gave  an  opportunity  to  make  some  interesting  comparisons,  elaborated  upon  in  this  study. 

Section  7.3  summarizes  observations  from  22  September  1993,  when  the  US  Department 
of  Energy  detonated  a  one  kiloton  conventional  explosion  at  the  Nevada  Test  Site  (NTS), 
in  an  experiment  named  the  Non-Proliferation  Experiment  (NPE).  The  experiment  was 
conducted  in  tire  context  of  future  CTBT/NPT  monitoring,  and  it  has  the  potential  of  pro¬ 
viding  data  useful  for  discrimination  between  chemical  and  nuclear  explosions.  The 
explosion  was  detected  automatically  at  three  of  the  northern/central  European  arrays: 
NORESS,  NORSAR  and  GERESS.  The  NPE  explosion  was  detonated  in  the  N-tunnel  in 
Rainier  Mesa  very  close  to  the  hypocenters  of  several  previous  nuclear  explosions,  some 
of  which  had  yields  similar  to  that  of  the  NPE  explosion.  In  particular,  we  compared  the 
NPE  to  the  nuclear  explosion  “Hunters  TVophy”  on  18  September  1992.  We  found  dial 
these  two  explosions  showed  very  similar  characteristics,  both  with  regard  to  signal  shape 
and  magnitude  (4. 14  for  “Hunters  Trophy”,  4. 10  for  NPE,  based  on  IMS  results).  With  an 
mb  of  4.4  as  given  by  USGS  to  “Hunters  Trophy”,  it  is  clear  that  NORSAR ’s  recordings 
of  the  NPE  event  confirm  that  this  one  kiloton  conventional  explosion  produced  ampli¬ 
tudes  somewhat  larger  than  those  expected  from  a  one  kiloton  fully  contained  nuclear 
explosion. 

Section  7.4  describes  a  generic  algorithm  for  accurate  determination  of  P-phase  arrival 
times.  This  is  a  continuation  of  our  studies  reported  mi  in  the  previous  NORSAR  Semian¬ 
nual  Technical  Summary.  In  that  report,  we  focused  upon  an  optimum  processing 
approach  for  one  particular  mining  area  (Khibiny  Massif).  We  have  now  developed  a 
generic  procedure  to  reestimate  the  onsets  of  all  types  of  first- arriving  P-phases,  assuming 
that  the  phase  identification  is  known.  By  applying  the  autoregressive  likelihood  tech¬ 
nique,  we  have  obtained  automatic  onset  times  of  a  quality  such  that  70%  of  the  automatic 
picks  are  within  0.1  s  of  the  best  manual  pick.  For  tire  onset  time  procedure  currently  used 
by  IMS,  die  corresponding  number  is  28%.  This  confirms  that  automatic  reestimation  of 
first-arriving  P-onsets  using  die  autoregressive  likelihood  technique  has  the  potential  of 
significantly  reducing  the  retiming  efforts  of  the  analyst. 
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Section  7  5  likewise  represents  a  follow-up  on  the  Intelligent  Post-Processing  (IPP)  of 
seismic  events  proposed  in  die  previous  Semiannual  Technical  Summary.  This  procedure 
has  been  shown  to  give  a  substantial  improvement  in  location  accuracy  when  applied  to 
seismic  events  in  the  Khibiny  Massif,  Kola  Peninsula.  In  this  paper,  we  compare  die  per¬ 
formance  of  the  analyst  using  the  Analyst  Review  Station  (ARS)  to  these  results.  As  part 
of  the  study,  we  estimate  the  uncertainty  in  analyst  time  picks  for  phases  from  various 
regional  arrays,  and  we  discuss  the  implication  of  these  uncertainties  in  terms  of  the  result¬ 
ing  effect  on  location  accuracy.  An  important  conclusion  inferred  from  this  work  is  that,  in 
many  cases,  location  accuracy  does  not  improve  when  ing  new  phase  readings  and 
applying  current  location  programs. 

Section  7.6  describes  an  application  of  Generalized  B camforming  (GBF)  applied  to  die  W. 
Caucasus  aftershock  sequence  of  29  April  1991.  This  sequence  occurred  during  the 
GSETT-2  main  phase,  and  the  GBP  results  could  therefore  be  compared  to  the  results  of 
the  four  Experimental  International  Data  Centers  (after  reprocessing).  The  GBF  associa¬ 
tion  process  reported  more  events  than  any  of  the  four  EIDCs,  and  also  had  the  most 
events  corresponding  to  the  reference  catalogue  (Starovoit  et  al,  1991).  In  addition,  the 
GBF  method  produced  17  reports  that  did  not  correspond  to  entries  in  Starovoit  et  al’s  bul¬ 
letin.  Each  of  the  EIDCs  also  had  events  in  this  category,  but  not  as  many  as  the  GBF  pro¬ 
cess.  One  event  reported  by  one  EIDC  and  confirmed  by  Starovoit  et  al’s  bulletin  was  not 
reported  by  the  GBF  method.  The  reason  was  that  the  event  had  only  two  valid  phases, 
and  thus  did  not  satisfy  our  GBF  detection  criterion.  On  the  other  hand,  the  GBF  reported 
4  confirmed  events  that  were  not  in  any  of  the  EIDC  bulletins.  The  study  concludes  that 
the  GBF  technique  provides  a  simple  and  rapid  way  to  associate  large  numbers  of  phases 
from  an  aftershock  sequence  with  a  very  low  false  alarm  rate.  In  fact,  the  GBF  aftershock 
processing  of  24  hours  of  data  for  the  day  in  question  (29  April  1991)  took  only  S  minutes 
on  a  SUN  sparcstation  2. 
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1  Summary 

This  Semiannual  Technical  Summary  describes  the  operation,  maintenance  and  research 
activities  at  the  Norwegian  Seismic  Array  (NORSAR),  the  Norwegian  Regional  Seismic 
Array  (NORESS),  the  Arctic  Regional  Seismic  Array  (ARCESS)  and  the  experimental 
Spitsbergen  regional  array  for  the  period  1  April  •  30  September  1993.  Statistics  are  also 
presented  for  additional  seismic  stations,  which  through  cooperative  agreements  with 
institutions  in  the  host  countries  provide  continuous  data  to  the  NORSAR  Data  Processing 
Center  (NIOC).  These  stations  comprise  the  Finnish  Experimental  Seismic  Array 
(FINESA),  the  German  Experimental  Seismic  Array  (GERESS),  and  an  experimental 
regional  seismic  array  in  Apatity,  Russia. 

This  Semiannual  Report  also  presents  statistics  from  operation  of  the  Intelligent  Monitor¬ 
ing  System  (IMS).  The  IMS  has  been  operated  in  an  experimental  mode,  and  the  perfor¬ 
mance  has  been  very  satisfactory.  Since  October  1991,  a  new  version  of  the  IMS  that 
accepts  data  from  an  arbitrary  number  of  arrays  and  single  3-component  stations  has  been 
operated. 

The  NORSAR  Detection  Processing  system  has  been  operated  throughout  the  period  with 
an  average  uptime  of  97.3%  as  compared  to  97.8%  for  the  previous  reporting  period.  A 
total  of  2290  seismic  events  have  been  reported  in  the  NORSAR  monthly  seismic  bulletin. 
The  performance  of  the  continuous  alarm  system  and  the  automatic  bulletin  transfer  by 
telex  to  AFTAC  has  been  satisfactory.  The  system  for  direct  retrieval  of  NORSAR  wave- 
farm  data  through  an  X.23  connection  has  been  used  successfully  for  acquiring  such  data 
by  AFTAC.  Processing  of  requests  for  full  NORSAR  and  regional  array  data  on  magnetic 
tapes  has  progressed  according  to  established  schedules.  There  have  been  no  modifica¬ 
tions  made  to  the  NORSAR  data  acquisition  system. 

On-line  detection  processing  and  data  recording  at  the  NORSAR  Data  Processing  Center 
(NDPC)  of  NORESS,  ARCESS  and  GERESS  data  have  been  conducted  throughout  the 
period.  FINESA  data  were  processed  until  16  May,  when  field  work  for  the  array  refur¬ 
bishment  started.  Data  from  two  experimental  small-aperture  arrays  at  sites  in  Spitsbergen 
and  Apatity,  Kola  Peninsula,  have  been  recorded  and  processed  in  an  experimental  mode. 
Moodily  processing  statistics  for  the  arrays  as  well  as  results  of  the  IMS  analysis  for  die 
reporting  period  are  given. 

Maintenance  activities  in  the  period  comprise  preventive/corrective  maintenance  in  con¬ 
nection  with  all  the  NORSAR  subarrays,  NORESS  and  ARCESS.  Other  activities  have 
involved  testing  of  die  NORSAR  communications  systems  and  work  in  connection  with 
the  experimental  small-aperture  arrays  in  Spitsbergen  and  Russia. 

Starting  1  October  1991,  an  effort  began  to  carry  out  a  complete  technical  refurbishment 
of  the  NORSAR  array.  This  project  is  funded  jointly  by  AFTAC,  ARPA  and  NFR.  During 
the  reporting  period,  we  have  continued  evaluation  and  laboratory  testing  of  technical 
options  for  field  instrumentation,  in  particular  state-of-the-art  A/D  converters,  data  acqui¬ 
sition  and  synchronization  devices.  During  the  reporting  period,  we  have  also  received  a 
complete  subarray  evaluation  unit  from  Science  Horizons,  and  this  unit  has  been  subjected 
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to  ft  long-term  operational  test  in  subarray  06C.  In  addition,  we  have  been  converting  the 
NORSAR  Detection/Event  Processing  software  to  be  compatible  with  our  UNIX-based 
workstations.  Our  detailed  status  regarding  this  work  is  reported  on  separately,  most 
recently  in  a  status  report  dated  16  August  1993. 

Summaries  of  six  scientific  contributions  are  presented  in  Chapter  7  of  this  report. 

Section  7.1  describes  an  experiment  in  continuous  threshold  monitoring  of  the  Lop  Nor, 
Chw,  nuclear  test  site.  The  monitoring  period  comprised  the  five  days  up  to  and  includ¬ 
ing  the  day  of  the  latest  nuclear  explosion  at  Lop  Nor  (5  October  1993).  NORESS, 
ARCESS  and  GERESS  regional  array  data  were  used.  We  have  been  compiling  daily  sta¬ 
tistics  of  all  peaks  on  the  threshold  diagram  exceeding  n^,  3.3  and  associating  these  peaks 
to  regional  or  teleseismic  events  whenever  possible.  Based  on  this  experiment,  it  appears 
that  this  threshold  (3.3),  which  is  about  0.73  mj,  units  higher  than  for  Novaya  Zemlya,  is 
representative  for  the  combined  threshold  monitoring  capabilities  (at  the  99  per  cent  level) 
of  these  three  regional  arrays. 

Section  7.2  presents  some  observations  for  the  =  3.9  Lop  Nor  explosion  on  S  October 
1993.  This  explosion  was  naturally  observed  on  the  threshold  plot  of  Section  7. 1  as  a  very 
sharp  peak.  The  contribution  summarizes  the  automatic  processing  results  for  all  six 
arrays  processed  at  Kjeller.  The  NORESS  array  has  the  best  signal-to-noise  ratio  (1376.2) 
for  this  event,  and  by  extrapolation  this  array  would  be  expected  to  have  a  detectable  sig¬ 
nal  for  an  event  from  Lop  Nor  about  2.3  magnitude  units  lower.  ARCESS  and  NORSAR 
also  show  outstanding  SNR.  The  velocity/azimuth  estimates  are  within  the  expected 
uncertainty  for  all  arrays.  It  is  also  noted  that  a  large  earthquake  and  aftershock  sequence 
occurred  in  the  Southern  Sinkiang  province  only  three  days  before  the  explosion,  and  this 
gave  an  opportunity  to  make  some  interesting  comparisons,  elaborated  upon  in  this  study. 

Section  7.3  summarizes  observations  from  22  September  1993,  when  the  US  Department 
of  Energy  detonated  a  one  kiloton  conventional  explosion  at  the  Nevada  Test  Site  (NTS), 
in  an  experiment  named  the  Non-Proliferation  Experiment  (NPE).  The  experiment  was 
conducted  in  the  context  of  future  CTBT/NPT  monitoring,  and  it  has  the  potential  of  pro¬ 
viding  data  useful  for  discrimination  between  chemical  and  nuclear  explosions.  The 
explosion  was  detected  automatically  at  three  of  the  northern/central  European  arrays: 
NORESS,  NORSAR  and  GERESS.  The  NPE  explosion  was  detonated  in  the  N-tunnel  in 
Rainier  Mesa  very  close  to  the  hypocenters  of  several  previous  nuclear  explosions,  some 
of  which  had  yields  similar  to  that  of  the  NPE  explosion.  In  particular,  we  compared  the 
NPE  to  the  nuclear  explosion  “Hunters  TVophy”  on  18  September  1992.  We  found  that 
these  two  explosions  showed  very  similar  characteristics,  both  with  regard  to  signal  shape 
and  magnitude  (4.14  for  “Hunters  TVophy”,  4. 10  for  NPE,  based  on  IMS  results).  With  an 
mj,  of  4.4  as  given  by  USGS  for  “Hunters  TVophy”,  it  is  clear  that  NORSAR ’s  recordings 
of  the  NPE  event  confirm  that  this  one  kiloton  conventional  explosion  produced  ampli¬ 
tudes  somewhat  larger  than  those  expected  from  a  one  kiloton  fully  contained  nuclear 
explosion. 

Section  7.4  describes  a  generic  algorithm  for  accurate  determination  of  P-phase  arrival 
times.  This  is  a  continuation  of  our  studies  reported  on  in  the  previous  NORSAR  Sernian- 
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nual  Technical  Summary.  In  that  report,  we  focused  upon  an  optimum  processing 
approach  for  one  particular  mining  area  (Khibiny  Massif).  We  have  now  developed  a 
generic  procedure  to  reestimate  the  onsets  of  all  types  of  fit  4-arriving  P-phases,  assuming 
that  the  phase  identification  is  known.  By  applying  the  autoregressive  likelihood  tech¬ 
nique,  we  have  obtained  automatic  onset  times  of  a  quality  such  that  70%  of  die  automatic 
picks  are  within  0.1  s  of  the  best  manual  pick.  Par  the  onset  time  procedure  currently  used 
by  IMS,  the  corresponding  number  is  28%.  This  confirms  that  automatic  reestimation  of 
first-arriving  P- onsets  using  the  autoregressive  likelihood  technique  has  the  potential  of 
significantly  reducing  the  retiming  efforts  of  the  analyst. 

Section  7.5  likewise  represents  a  follow-up  on  the  Intelligent  Post-Processing  (IPP)  of 
seismic  events  proposed  in  the  previous  Semiannual  Technical  Summary.  This  procedure 
has  been  shown  to  give  a  substantial  improvement  in  location  accuracy  when  applied  to 
seismic  events  in  the  Khibiny  Massif,  Kola  Peninsula.  In  this  paper,  we  compare  the  per¬ 
formance  of  the  analyst  using  the  Analyst  Review  Station  (ARS)  to  these  results.  As  part 
of  the  study,  we  estimate  the  uncertainty  in  analyst  time  picks  for  phases  from  various 
regional  arrays,  and  we  discuss  the  implication  of  these  uncertainties  in  terms  of  the  result¬ 
ing  effect  on  location  accuracy.  An  important  conclusion  inferred  from  this  work  is  that,  in 
many  cases,  location  accuracy  does  not  improve  when  adding  new  phase  readings  and 
applying  current  location  programs. 

Section  7.6  describes  an  application  of  Generalized  Beamforming  (GBF)  applied  to  the  W. 
Caucasus  aftershock  sequence  of  29  April  1991.  This  sequence  occurred  during  the 
GSETT-2  main  phase,  and  the  GBF  results  could  therefore  be  compared  to  the  results  of 
the  four  Experimental  International  Data  Centers  (after  reprocessing).  The  GBF  associa¬ 
tion  process  reported  more  events  than  any  of  the  four  EIDCs,  and  also  had  the  most 
events  corresponding  to  the  reference  catalogue  (Starovoit  et  al,  1991).  In  addition,  the 
GBF  method  produced  17  reports  that  did  not  correspond  to  entries  in  Starovoit  et  al’s  bul¬ 
letin.  Each  of  the  EIDCs  also  had  events  in  this  category,  but  not  as  many  as  the  GBF  pro¬ 
cess.  One  event  reported  by  one  EIDC  and  confirmed  by  Starovoit  et  al’s  bulletin  was  not 
reported  by  the  GBF  method.  The  reason  was  that  the  event  had  only  two  valid  phases, 
and  thus  did  not  satisfy  our  GBF  detection  criterion.  On  the  other  hand,  the  GBF  reported 
4  confirmed  events  that  were  not  in  any  of  the  EIDC  bulletins.  The  study  concludes  thai 
the  GBF  technique  provides  a  simple  and  rapid  way  to  associate  large  numbers  of  phases 
from  an  aftershock  sequence  with  a  very  low  false  alarm  rate.  In  fact,  the  GBF  aftershock 
processing  of  24  hours  of  data  for  the  day  in  question  (29  April  1991)  took  only  5  minutes 
on  a  SUN  sparcstation  2. 
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2  NORSAR  Operation 

2.1  Detection  Processor  (DP)  operation 

There  have  been  82  breaks  in  the  otherwise  continuous  operation  of  the  NORSAR  online 
system  within  the  current  6-month  reporting  interval.  The  uptime  percentage  for  the 
period  is  97.3  as  compared  to  97.8  for  the  previous  period. 

Fig.  2.1.1  and  the  accompanying  Table  2.1.1  both  show  the  daily  DP  downtime  for  the 
days  between  1  April  and  30  September  1993.  The  monthly  recording  times  and  percent¬ 
ages  are  given  in  Table  2.1.2. 


The  breaks  can  be  grouped  as  follows: 


a) 

Hardware  failure 

39 

b) 

Stops  related  to  program  work  or  error 

0 

c) 

Hardware  maintenance  stops 

5 

d) 

Power  jumps  and  breaks 

0 

e) 

TOD  error  correction 

12 

f) 

Communication  lines 

26 

The  total  downtime  for  the  period  was  127  hours  and  27  minutes.  The  mean-time- 
be tween-failures  (MTBF)  was  2.8  days,  as  compared  to  3.6  for  the  previous  period. 


J.  Torstveit 
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Fig.  2.1.1.  Detection  Processor  uptime  for  April  (top),  May  (middle)  and  June  (bottom) 
1993. 
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Fig.  2.1.1,  Detection  Processor  uptime  for  July  (top),  August  (middle)  and  September 
(bottom)  1993. 
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Date 

Time 

21 

Jun  2045  - 

22  Jun 

-  0032 

30  Jun 

0700  -0744 

03  Jul 

1150  -  1427 

06  Jul 

0304  -  0437 

06  Jul 

0728  - 

07  Jul 

-  1514 

08  Jul 

1500  -  1540 

08  Jul 

2200  -  2306 

09  Jul 

0532  -  0613 

09  Jul 

1104  -  1143 

09  Jul 

1441  -  1502 

09  Jul 

2200  -  2237 

10  Jul 

2200  -  2239 

11  Jul 

1407  -  1442 

11  Jul 

2200  - 

12  Jul 

-  0531 

12  Jul 

0744  -  1307 

10  Aug 

0821  -  1514 

29  Aug 

1000  - 

30  Aug 

-  1211 

08  Sep 

1518  - 

09  Sep 

-  0428 

10  Sep 

0821  -  1030 

18  Sep 

0817  -  1024 

20  Sep 

1208  -  1316 

22  Sep 

0804  -  0825 

22  Sep 

2230  -  2317 

27  Sep 

1719  -  1759 

Table  2.1.1.  The  major  downtimes 

Cause 

Line  failure 
Line  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  failure 
Hardware  service 
Line  failure 
Line  failure 
Line  failure 
Line  failure 
Line  failure 
Hardware  failure 
Hardware  failure 
Hardware  service 
Hardware  failure 
Hardware  failure 
Hardware  failure 


the  period  1  April  -  30  September  1993. 
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Month 

DP  Uptime 
Hours 

DP  Uptime 
% 

No.  of 

DP  Breaks 

No.  of 
Days  with 
Breaks 

DPMTBF* 

(days) 

Apr  93 

719.29 

99.93 

4 

4 

4.6 

May  93 

743.23 

99.92 

6 

4 

3.4 

Jun  93 

715.29 

99.37 

2 

3 

7.9 

Jul93 

686.46 

92.40 

37 

14 

0.6 

Aug  93 

701.15 

95.47 

7 

8 

2.8 

Sep  93 

698.11 

96.97 

26 

18 

0.8 

97.34 

82 

_ 

53 

2.8 

•Mean-time-between-failures  =  total  uptime/oo.  of  up  intervals. 


Table  2.12.  Online  system  performance,  1  April  -  30  September  1993. 
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|  2.2  Array  communications 

General 

Table  2.2.1  reflects  the  performance  of  the  communications  system  throughout  the  report¬ 
ing  period. 

The  following  lists  the  most  prominent  events  which  have  affected  individual  and/or  all 
systems  simultaneously: 

•  Loss  of  synch 

•  Damaged  cables 

•  Scheduled  communications  cable  work  (NTA) 

•  Unspecified  wait  related  to  communications  line  (NTA) 

•  Line  unit  /CTV  modem 

•  Timing  problems  after  installation  of  digital  central 

,  •  Failing  DACCS  (Digital  Access  Cross  Connection  System/NTA) 

•  Scheduled  maintenance  DACCS  Lilies  trom/Hamar 

•  Lightning  in  CTV  area  (spec.  02B,  04C,  06C) 

•  Modcomp  memory  chip  failure 

•  IBM  disk  3656  controller  failure  (loss  of  power) 

•  2701  IBM  Data  Adapter  failure 

•  2701 -Modcomp  data  transfer  stop 

Detailed  Summary 

April  (weeks  13-17),  29.3-2.5.93 

The  subaarray  02C  was  turned  off  26  April  at  1349  hrs  due  to  bad  status  in  the  communi¬ 
cations  system  (NO  DATA).  A  modcomp  restart  27  April  reinitiated  the  operation. 

A  NORSAR  stop  occurred  29  April  at  0903  hrs,  probably  due  to  an  irregularity  between 
4  the  2701  adapter  and  modcomp.  The  system  was  restarted  at  0922  hrs. 

Average  outages  for  April,  individual  weeks 


Week  13  (all) 

0.004% 

Week  14  (all) 

0.014% 

Week  15  (all) 

0.0003% 

Week  16  (all) 

0.005% 

Week  17  (-02C) 

0.0009% 

i 


May  (weeks  18-21),  3-30.5.93 

In  the  period  the  three  subarrays  02B,  02C  and  06C  have  been  affected.  Modcomp  restarts 
reinitiated  the  systems  again. 
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Average  outages  in  May,  individual  weeks: 

Week  18  (all)  0.004% 

Week  19  (all)  0.036% 

Week  20  (all)  0.005% 

Week  21(all)  0.117% 

June  (weeks  22-26),  31.5-4.7.93 

All  systems  were  affected  in  the  period.  04C  was  down  1  June  between  2000  and  2242  hrs 
GMT  in  connection  with  scheduled  cable  work. 

On  5  June  data  was  lost  from  02B  kl  1701  GMT  caused  by  lightning.  Electricians  from 
Rena  Power  Plant  replaced  fuses  1240  hrs  GMT  7  June. 

02B  was  affected  again  on  18  June  between  0814  and  1230  hrs  GMT.  NTA  carried  out 
cable  rerouting. 

On  21  June  02C  data  disappeared  after  NTA/Hamar  had  replaced  the  original  analog  cen¬ 
tral  at  the  NTA  Sjursoen  premises  with  a  digital  central  (PCM).  In  addition,  the  operation 
also  involved  a  change  in  the  routing  of  data. 

The  same  day  NTA/Lillestrem  carried  out  scheduled  maintenance  of  a  DACCS  in  the 
NTA/Lillestrom  premises,  affecting  all  the  NORSAR  subarrays  and  the  64  Kbit  Spits¬ 
bergen  line.  All  the  subarrays  were  masked  2045  hrs  GMT. 

On  22  June,  at  0038  hrs  GMT,  all  the  subarrays  (-01  A,  02C)  were  demasked. 

On  24  June  01A  resumed  operation  after  replacing  a  line  unit  in  the  CTV  modem. 

On  30  June  the  communications  between  the  2701  Data  Adapter  and  modem/modcomp 
stopped  for  approximately  44  minutes. 

NORSAR  time  was  set  back  1  second  on  1  July. 

On  2  July  a  modcomp  memory  chip  failed  and  was  replaced.  Downtime  was  approxi¬ 
mately  2  hrs  35  minutes. 

Average  outage  in  June,  individual  weeks: 


Week  22  (-02B)  0.002% 

Week  23  (-02B)  0.002% 

Week  24  (all)  0.002% 

Week  25  (-01A.02C)  0.003% 

Week  26  (-02C)  1.488% 


July  (weeks  27-30),  5.7-1.8.93 

3  July  the  modcomp  failed.  A  memory  chip  was  replaced  and  the  system  was  started 
again.  Approximate  downtime  was  2.5  hrs. 
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On  6  July  an  IBM  disk  controller  (channel  1)  lost  power  approximately  0304  hrs.  Started 
again  0437  hrs. 

The  2701  Data  Adapter  also  failed  on  6  July.  After  the  +3  Volt  power  supply  was  repaired, 
the  system  was  started  again.  NORSAR  downtime  was  28  hrs. 

On  12  July  Mem  ore  x  replaced  the  +5  Volt  power  supply  in  the  disk  controller  3656  (chan¬ 
nel  1).  Between  7  and  12  July  the  controller  stopped  28  times. 

19*20  July  the  04C  communications  line  was  down  approximately  16.5  hrs. 

At  the  end  of  July  02C  was  still  down 

Average  outages  in  July,  individual  weeks: 


Week  27  (-02C)  0.0005% 

Week  28  (-02C)  0.0009% 

Week  29  (-02C.04C)  0.0008% 

Week  30  (-02C.04Q  0.002% 


August  (weeks  31-34),  2-29.8.93 

The  02C  communications  problems  started  21  June  when  the  existing  communications 
equipment  was  replaced  (modified)  with  PCM  equipment  in  the  Sjusjeen  area,  thereby 
introducing  timing  problems.  The  first  step  in  trying  to  solve  the  problem  was  to  replace 
the  NDPC/CTV  modems  as  quickly  as  possible. 

02B  was  down  between  4  and  12  August  due  to  a  broken  cable  between  Kjeller  and 
Lillestrem. 

06C  was  down  between  7  and  9  August,  probably  caused  by  synch  problems. 

On  10  August  all  the  subarrays  were  affected  by  a  failing  DACCS  in  the  NTA  premises/ 
Hamar. 

On  19  August  06C  was  again  affected  and  resumed  operation  after  a  modcomp  restart  23 
August. 

On  26  August  NTA/Hamar  carried  out  unspecified  work  related  to  the  01 A  communica¬ 
tions  system.  Two  days  later,  28  August,  the  performance  of  the  01A  communications  sys¬ 
tem  gradually  was  reduced.  On  30  August  01 A  was  declared  down  and  masked. 

Average  outages  in  August,  individual  weeks: 


Week  31  (-02B,02C,06Q  0.001% 

Week  32  (-02Q  3.610% 

Week  33  (-02C.06C)  0.003% 

Week  34  (-02Q  0.050% 
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September  (weeks  35-39),  30.8-3.10.93 

01  A,  which  went  down  28  August  due  to  a  bad  cable,  was  restarted  13  September  after 
repair. 

Between  8  September  1306  hrs  and  9  September  0428  hrs  all  communications  systems 
were  affected  (13.5  hrs). 

Reduced  01B  performance  started  10  September  when  stop  messages  were  received  per¬ 
taining  to  the  communications  lines. 

7  and  8  September  NTA/Lillestram  and  Hamar  tried  to  solve  the  02C  communications 
problems.  The  problem  started  when  a  new  PCM-system  was  introduced  21  June  93. 
Their  first  initiative  was  to  install  new  modems  hoe  at  Kjeller  and  in  die  CTV.  The  test 
revealed  error-free  transmission  in  both  directions  modem  to  modem,  but  failed  when  the 
remaining  NORS  AR  equipment  was  attached.  The  conclusion  was  tuning  problems.  New 
initiatives  will  be  taken. 

The  01B  communications  cable  was  repaired,  according  to  NTA/Hamar,  on  13  September. 
An  attempt  to  start  the  line  failed. 

Between  9  and  20  September  there  were  a  large  number  of  modcomp  restarts  in  connec¬ 
tion  with  the  misalignment  of  data  related  to  clock  pulses. 

On  22  September  two  stops  occurred,  in  both  cases  because  the  data  transfer  between  the 
2701  and  the  high-speed  modems  had  stopped. 

Also  on  23  and  29  September  and  1  October  there  were  misalignments  of  data  related  to 
clock  pulses.  In  order  to  resynch  die  06C,  the  modcomp  was  restarted  28  September. 

There  seems  to  be  a  connection  between  misalignment  of  data  related  to  clock  pulses  and 
the  DACCS  installed  in  die  NTA  premises  in  Lillestrpm  and  Hamar.  We  have  observed 
data  misalignment  both  before  and  after  an  announced  DACCS  corrective  maintenance  30 
September  (which  failed).  Finally,  we  lost  data  between  0915  hrs  and  1130  hrs  in  connec¬ 
tion  with  DACCS  replacment  at  NTA/Lillestram. 

Average  outages  in  September,  individual  weeks: 


Week  35  (-01A,  02C)  0.003% 

Week  36  (-01  A,01B,02Q  0.56% 

Week  37  (N/A) 

Week  38  (-01B.02C)  0.005% 

Week  39  (-01B,02C,06C)  1 .49  % 


OA.  Hansen 
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Snbsurays 

Apr  (5) 

May  (4) 

Jun  (5) 

Jul  (4) 

Aug  (4) 

Sep  (5) 

Average 

29.3-2.5.93 

3-303.93 

31.5-4.7.93 

5.7-1.8.93 

2-29.8.93 

30.8- 

3.10.93 

1/2  year 

01A 

0.005 

0.002 

ilSQfli 

0.0007 

0.024u) 

N/A 

01B 

0.004 

0.0006 

0.0006® 

0.001 

0.001 12) 

N/A 

0.000724) 

EB 

0.009 

0.009 

0.0004S) 

0.001 

0.01313) 

0.0121® 

0.007 

02C 

0.020l) 

0.012 

0.006® 

N/A 

N/A 

N/A 

0.01 225) 

03C 

0.005 

0.0006 

o.warf1 

0.0004 

0.000414> 

0.00417) 

0.001 

04C 

0.001 

0.002 

0.001® 

0.001 10) 

0.0021® 

0.0021® 

0.001 

06C 

0.001 

03402) 

0.0019) 

0.001 

N/A 

0.00619) 

0.0102® 

AVER 

0.005 

0.052 

0.001 

O-OOOS20) 

0.008?' * 

0.00622' 

0.005 

Figures  representing  error  rate  (in  per  cent)  followed  by  number  1).  2).  etc.,  tre  related  to  legend  below. 


Table  2.2.1.  Communications  performance.  The  numbers  represent  error  rates  in  per  cent 
based  on  total  transmitted  frames/week  (29  Mar  -  3  Oct  93). 


5).  1®.  13) 

2).  3).  ®.ll).  13) 
14).  16).  18).  19) 
1).4).7).8).9).  17) 
25) 

23).  24) 


26) 

20) 

21) 

22) 


avenge  2  weeks  (2425/2728) 

avenge  3  weeks  (192021/222324/313334/3526.37/352738) 

avenge  4  weeks  (13.14.15.16/22232425/25262728) 
avenge  3  months  (Apr.  May,  Jun) 
avenge  5  months  (Apr;  May,  Jun,  Jul,  Aug) 
avenge  5  months  (Apr,  May.  Jun,  Jul  Sep) 
avenge  6  subarrays  (01A.02B,  03C-060 
average  5  subways  (01A-02B,  03C,  040 
average  4  subamys  (02B,  03C-06Q 
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23  NORSAR  Event  Detection  operation 

In  Txblc  2.3.1  some  monthly  statistics  of  the  Detection  and  Event  Processor  operation  are 
given.  The  table  lists  the  total  number  of  detections  (DPX)  triggered  by  the  on-line  detec¬ 
tor,  the  total  number  of  detections  processed  by  the  automatic  event  processes'  (EPX)  and 
tile  total  number  of  events  accepted  after  analyst  review  (teleseismic  phases,  core  phases 
and  total). 


Total 

DPX 

Total 

EPX 

Accepted  events 

Sum 

Daily 

P-phase8 

Core 

Phases 

Apr  93 

10725 

1398 

248 

76 

324 

10.8 

May  93 

4525 

830 

305 

92 

397 

12.8 

Jun93 

7825 

1214 

292 

130 

422 

14.1 

Jul93 

6750 

1146 

342 

68 

410 

13.2 

Aug  93 

7900 

1232 

285 

54 

339 

10.9 

Sep  93 

9575 

1418 

285 

113 

398 

13.3 

1757 

533 

2290 

12.5 

Table  2.3. 1 .  Detection  and  Event  Processor  statistics,  1  Apr  -  30  Sep  1993. 


NORSAR  Detections 

The  number  of  detections  (phases)  reported  by  the  NORSAR  detector  during  day  091 
through  day  273, 1993,  waa  47,524,  giving  an  average  of  263  detections  per  processed  day 
(183  days  processed).  Table  2.3.2  shows  daily  and  hourly  distribution  of  detections  for 
NORSAR. 
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Table  23.2.  Daily  and  hourly  distribution  of  NORS  AR  detections.  For  each  day  is  shown  number  of  detections  within  each 
hour  of  the  day  and  number  of  detections  for  that  day.  The  end  statistics  give  total  number  of  detections  distributed  for  each 
hour  and  the  total  sum  of  detections  during  the  period.  The  averages  show  number  of  processed  days,  hourly  distribution  and 
average  per  processed  day.  (Page  4  of  4) 
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3  Operation  of  regional  arrays 

3.1  Recording  of  NORESS  data  at  NDPC,  Kjeller 

Table  3.1.1  lists  the  main  outage  times  and  reasons. 

The  average  recording  time  was  99.94%  as  compared  to  99.96%  during  the  previous 
reportin,,  period. 


Date 

Time 

Cause 

14  Jul 

2245  -  2254 

Transmission  line  failure 

03  Aug 

1054  -  1224 

Hardware  maintenance 

26  Sep 

0100  -  0200 

Software  failure 

Table  3.1.1.  Interruptions  in  recording  of  NORESS  data  at  NDPC,  1  April  -  30  September 
1993. 


Monthly  uptimes  for  the  NORESS  on-line  data  recording  task,  taking  into  account  all  fac¬ 
tors  (field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were 
as  follows: 


April 

100.00 

May 

100.00 

June 

100.00 

July 

99.97 

August 

99.79 

September 

99.86 

Fig.  3.1.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
NORESS  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period. 

J.  Torstveit 
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Fig.  3.1.1,  NORESS  data  recording  uptime  for  Apirl  (top).  May  (middle)  and  June  (bot¬ 
tom)  1993. 
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Fig.  3.1.1.  (cont.)  NORESS  data  recording  uptime  for  July  (top),  August  (middle)  and 
September  (bottom)  1993. 
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3.2  Recording  of  ARCESS  data  at  NDPC,  Kjelier 

Table  3.2.1  lists  the  main  outage  times  and  reasons. 

The  average  recording  time  was  99.19%  as  compared  to  99.61%  for  the  previous  reporting 


period. 

Date 

Time 

Cause 

09  Apr 

1953 

-  2054 

Software  failure 

26  Apr 

0231 

-  0344 

Satellite  link  failure 

10  May 

1028 

-  1300 

Powerline  work  HUB 

11  May 

1039 

-  1235 

Powerline  work  HUB 

12  May 

1031 

-  1244 

Powerline  work  HUB 

13  May 

1041 

-  1227 

Powerline  work  HUB 

04  Jun 

2045 

Hardware  failure  HUB 

05  Jun 

-  0706 

Hardware  failure  HUB 

05  Jun 

1035 

-  1108 

Hardware  failure  NDPC 

06  Jul 

1049 

-  1348 

Hardware  failure  NDPC 

09  Jul 

0107 

-  0515 

Satellite  link  failure 

13  Jul 

0739 

-  0800 

Hardware  failure  NDPC 

10  Aug 

0707 

-  1038 

Satellite  link  service 

Table  3.2.1 .  The  main  interruptions  in  recording  of  ARCESS  data  at  NDPC,  1  April  -  30 
Septemer  1993. 

Monthly  uptimes  for  the  ARCESS  on-line  data  recording  task,  taking  into  account  all  fac¬ 
tors  (field  installations,  transmissions  line,  data  center  operation)  affecting  this  task  were 
as  follows: 


April 

99.68% 

May 

98.82% 

June 

98.23% 

July 

98.98% 

August 

99.48% 

March 

99.96% 

Pig.  3.2.1.  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
ARCESS  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period. 

J.  Torstveit 
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Fig.  &2.I.  ARCESS  data  recording  uptime  for  July  (top),  August  (middle)  and  September 
(bottom)  1993. 
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33  Recording  of  FINESA  data  at  NDPC,  Kjeller 


|  The  HUB  field  system  was  hit  by  lightning  cm  May  16  and  not  brought  back  into  operation 

|  again,  as  work  related  to  a  planned  refurbishment  started  shortly  afterwards. 

i 

] 


Date 

Time 

Cause 

11  Apr 

1056  - 

HUB  failure 

12  Apr 

-  0111 

HUB  failure 

16  Apr 

0832  -  1027 

Thmsmission  line  test 

08  May 

1404  -  1508 

Transmission  line  failure 

08  May 

1524  -  1544 

Transmission  line  failure 

08  May 

1551  -  1610 

Transmission  line  failure 

16  May 

1342  - 

HUB  failure,  field  system  upgrade 

Htble  3.3.1.  The  main  interruptions  in  recording  of  FINESA  data  at  NDPC,  1  April  -  16 
May  1993. 

Monthly  uptimes  for  the  FINESA  on-line  data  recording  task,  taking  into  account  all  fac¬ 
tors  (field  installations,  transmission  lines,  data  center  operation)  affecting  this  task  were 
as  follows: 


April 

97.68% 

May 

49.99% 

June 

0.00% 

July 

0.00% 

August 

0.00% 

September 

0.00% 

Fig.  3.3.1  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
FINESA  data  in  our  tape  archive,  on  a  day-by-day  basis,  for  the  reporting  period. 

J.  Torstveit 
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3.4  Recording  of  Spitsbergen  data  at  NDPC,  Kjeller 

The  average  recording  time  was  99.16%. 

The  main  reasons  for  downtime  follow: 


Date 

Time 

12  May 

1955  - 

2036 

21  May 

1944  - 

2026 

21  May 

2131  - 

2313 

23  May 

2223  - 

2300 

21  Jun 

2133  - 

2316 

30  Jun 

0743  - 

0827 

01  Jul 

1313  - 

1427 

01  Jul 

1518  - 

1612 

06  Jul 

1750  - 

1853 

06  Jul 

1909  - 

2002 

30  Jul 

1900  - 

2012 

02  Aug 

2000  - 

2141 

10  Aug 

1146  - 

1242 

10  Aug 

1939  - 

2010 

11  Aug 

0919  - 

1023 

12  Aug 

1555  - 

1647 

03  Sep 

0852  - 

1136 

06  Sep 

0842  - 

1000 

22  Sep 

0814  - 

0938 

24  Sep 

0439  - 

0600 

Cause 

Communication  line  failure 
Communication  line  failure 
Communication  line  failure 
Communication  line  failure 
Communication  line  failure 
Communication  line  failure 
Communication  line  failure 
Communication  line  failure 
Communication  line  failure 
Communication  line  failure 
Communication  line  failure 
Communication  line  failure 
Communication  line  failure 
Communication  line  failure 
Communication  line  failure 
Communication  line  failure 
Maintenance  work  at  array  site 
Maintenance  work  at  array  site 
Communication  line  failure 
Communication  line  failure 


Monthly  uptimes  for  the  Spitsbergen  online  data  recording  task,  taking  into  account  all 
factors  (field  installations,  transmission  line,  data  center  operation)  affecting  this  task  were 
as  follows: 


April 

99.91% 

May 

99.07% 

June 

99.43% 

July 

9858% 

August 

98.97% 

September 

99.00% 

Fig.  3.4. 1  shows  the  uptime  for  the  data  recording  task,  or  equivalently,  the  availability  of 
Spitsbergen  data  in  our  tape  archive,  on  a  day-by-day  basis  for  the  reporting  period. 

J.Tbrstveit 
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33  Event  detection  operation 

This  section  reports  results  from  one-array  automatic  processing  using  signal  processing 
recipes  and  “ronapp”  recipes  for  the  ep  program  (NORSAR  Sd.  Rep.  No  2-88419). 

Three  systems  are  in  parallel  operation  to  associate  detected  phases  and  locate  events: 

1.  The  ep  program  with  “ronapp”  recipes  is  operated  independently  on  each  array  to 
obtain  simple  one-array  automatic  solutions. 

2.  The  Generalized  Beamfcrming  method  (GBF)  (see  Ringdal  and  Kvtema  ( 1989),  “A 
mulitchannel  processing  approach  to  real  time  network  detection,  phase  association 
and  threshold  monitoring”,  Bull.  Seism.  Soc.  Am.  79, 6, 1927-1940)  processes  the 
regional  arrays  jointly  and  presents  1  oca  turns  of  regional  events. 

3.  The  IMS  system  is  operated  on  the  same  set  of  arrivals  as  ep  and  GBP  and  reports 
also  teleseismic  events  in  addition  to  regional  ones. 

IMS  results  are  reported  in  sections  3.6  ami  3.7  and  GBP  results  in  section  3.8. 

In  addition  to  these  three  event  association  processes,  we  are  running  test  versions  of  the 
so-called  Threshold  Monitoring  (TM)  process.  This  is  a  process  that  monitors  the  seismic 
amplitude  level  at  the  regional  arrays  continuously  in  time  to  estimate  the  upper  magni¬ 
tude  limit  of  an  event  that  might  go  undetected  by  the  network.  The  current  TM  process  is 
beamed  to  several  sites  of  interest,  including  the  Novaya  Zemlya  test  site.  Simple  displays 
of  so-called  threshold  curves  reveal  instants  of  particular  interest;  i.e.,  instants  when 
events  above  a  certain  magnitude  threshold  may  have  occurred  in  the  target  region. 
Results  from  the  three  processes  described  above  are  used  to  help  resolve  what  actually 
happened  during  these  instances. 

NORESS  detections 

The  number  of  detections  (phases)  reported  from  day  091, 1993,  through  day  273, 1993, 
was  31,905,  giving  an  average  of  174  detections  per  processed  day  (183  days  processed). 

Table  3.5.1  shows  daily  and  hourly  distribution  of  detections  for  NORESS. 

Events  automatically  located  by  NORESS 

During  days  091, 1993,  through  273, 1993, 2113  local  and  regional  events  were  located  by 
NORESS,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  11.5  events  per  processed  day  (183  days  processed).  68%  of  these  events  are  within 
300  km,  and  89%  of  these  events  are  within  1000  km. 

ARCESS  detections 

The  number  of  detections  (phases)  reported  during  day  091, 1993,  through  day  273  1993, 
was  72,987,  giving  an  average  of  401  detections  per  processed  day  (182  days  processed). 
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Table  3.S.2  shows  daily  and  hourly  distribution  of  detections  for  ARCESS. 

Events  automatically  located  by  ARCESS 

During  days  091, 1993,  through  273, 1993, 3917  local  and  regional  events  were  located  by 
ARCESS,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
21.3  events  per  processed  day  (182  days  processed).  51%  of  these  events  are  within 
300  km,  and  88%  of  these  events  are  within  1000  km. 

FINES  A  detections 

The  number  of  detections  (phases)  reported  during  day  091, 1993,  through  day  136, 1993, 
was  23,403,  giving  an  average  of  S09  detections  per  processed  day  (46  days  processed). 

Table  3.5.3  shows  daily  and  hourly  distribution  of  detections  for  FINESA. 

Events  automatically  located  by  FINESA 

During  days  091, 1993,  through  136, 1993, 999  local  and  regional  events  were  located  by 
FINESA,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  21.7  events  per  processed  day  (46  days  processed)  52%  of  these  events  are  within 
300  km,  and  81%  of  these  events  are  within  1000  km. 

GERESS  detections 

The  number  of  detections  (phases)  reported  from  day  091, 1993,  through  day  273, 1993, 
was  34,424,  giving  an  average  of  190  detections  per  processed  day  (181  days  processed). 

Table  3.5.4  shows  daily  and  hourly  distribution  of  detections  for  GERESS . 

Events  automatically  located  by  GERESS 

During  days  091, 1993,  through  273, 1993, 3444  local  and  regional  events  were  located  by 
GERESS,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an  average 
of  19.0  events  per  processed  day  (181  days  processed).  76%  of  these  events  are  within 
300  km,  and  91%  of  these  events  are  within  1000  km. 


Apatity  array  detections 

The  number  of  detections  (phases)  reported  from  day  091, 1993,  through  day  273, 1993, 
was  159,561,  giving  an  average  of  891  detections  per  processed  day  (179  days  processed). 

As  described  in  earlier  reports,  die  data  from  the  Apatity  array  are  transferred  by  one-way 
(simplex)  radio  links  to  Apatity  city.  The  transmission  suffers  from  radio  disturbances  that 
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results  in  a  large  number  of  small  data  gaps  and  spikes  in  the  data.  Although  die  communi¬ 
cation  protocol  may  correct  such  errors  by  requesting  retransmission  of  data,  this  cannot 
presently  be  done  at  Apatity.  For  such  error  corrections,  a  two-way  radio  link  is  needed 
(duplex  radio).  However,  it  should  be  noted  that  noise  from  cultural  activities  and  from  the 
nearby  lakes  cause  most  of  die  unwanted  detections.  These  unwanted  detections  are  “fil¬ 
tered”  in  the  signal  processing,  as  they  give  seismic  velocities  that  are  outside  accepted 
limits  for  regional  and  teleseism  ic  phase  velocities. 

Table  3.5  .5  shows  daily  and  hourly  distribution  of  detections  for  the  Apatity  array. 


Events  automatically  located  by  the  Apatity  array 

During  days  091, 1993,  through  273, 1993, 2142  local  and  regional  events  were  located  by 
the  Apatity  array,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives  an 
average  of  12. 1  events  per  processed  day  (179  days  processed).  41%  of  these  events  are 
within  300  km,  and  72%  of  these  events  are  within  1000  km. 


Spitsbergen  array  detections 

The  number  of  detections  (phases)  reported  from  day  091, 1993,  through  day  273, 1993, 
was  58,034,  giving  an  average  of  317  detections  per  processed  day  (183  days  processed). 

Table  3.5.6  shows  daily  and  hourly  distribution  of  detections  for  the  Spitsbergen  array. 

Events  automatically  located  by  the  Spitsbergen  array 

During  days  091, 1993,  through  273,1993, 3880  local  and  regional  events  were  located  by 
the  Spitsbergen  array,  based  on  automatic  association  of  P-  and  S-type  arrivals.  This  gives 
an  average  of  23. 1  events  per  processed  day  (183  days  processed).  34%  of  these  events  are 
within  300  km,  and  72%  of  these  events  are  within  1000  km. 
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Table  3.5.2,  Daily  and  hourly  distribution  of  ARCESS  detections.  For  each  day  is  shown  number  of  detections  within 
each  hour  of  die  day,  and  number  of  detections  for  that  day.  The  end  statistics  give  total  number  of  detections 
distributed  for  each  hour  and  the  total  sum  of  detections  during  the  period.  The  averages  show  number  of  pro¬ 
cessed  days,  hourly  distribution  and  average  per  processed  day.  (Page  4  of  4) 
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3.6  IMS  operation 

The  Intelligent  Monitoring  System  (IMS)  was  installed  at  NORSAR  in  December  1989 
and  was  operated  at  NORS  AR  from  1  January  1990  for  automatic  processing  of  data  from 
ARCESS  and  NORESS.  A  second  version  of  IMS  that  accepts  data  from  an  arbitrary 
number  of  arrays  and  single  3-component  stations  was  installed  at  NORSAR  in  October 
1991,  and  regular  operation  of  the  system  comprising  analysis  erf  data  from  the  4  arrays 
ARCESS,  NORESS,  FINESA  and  GERESS  started  on  IS  October  1991.  As  apposed  to 
the  first  version  of  IMS,  the  one  in  current  operation  also  locates  events  at  teleseismic  dis¬ 
tance. 

On  14  December  1992,  phase  detections  from  the  Apatity  array  were  included  in  the  auto¬ 
matic  phase  association.  The  phase  detections  from  the  Spitsbergen  array  were  made 
available  to  die  analysts  on  5  February  1993.  These  detections  are  not  used  in  the  auto¬ 
matic  phase  association,  but  can  be  added  manually  during  analysis. 

The  operational  stability  of  IMS  has  been  very  good  during  die  reporting  period.  In  fact 
the  IMS  event  processor  (pipeline)  has  had  no  downtime  of  its  own;  i.e.,  all  data  available 
to  IMS  have  been  processed  by  IMS. 

Phase  and  event  statistics 

Table  3.6. 1  gives  a  sumary  of  phase  detections  and  events  declared  by  IMS .  From  top  to 
bottom  the  table  gives  the  total  number  of  detections  by  the  IMS,  the  number  of  detections 
that  are  associated  with  events  automatically  declared  by  the  IMS,  the  number  of  detec¬ 
tions  that  are  not  associated  with  any  events,  die  number  of  events  automatically  declared 
by  the  IMS,  die  total  number  of  events  defined  by  the  analyst,  and  finally  the  number  of 
events  accepted  by  the  analyst  without  any  changes  (i.e.,  from  the  set  of  events  automati¬ 
cally  declared  by  the  IMS) 

U.  Baadshaug 
B.Ferstad 
B.Kr.  Hokland 
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Apr  93 

May  93 

Jun93 

Jul93 

Aug  93 

Sep  93 

Total 

Phase  detections 

S4641 

78795 

67666 

57248 

59054 

57115 

374519 

8757 

7376 

6206 

6395 

6785 

6749 

42268 

-Unassodated  phases 

45884 

71419 

61460 

50853 

52269 

50366 

332251 

Events  automatically 
declared  by  IMS 

2500 

2206 

1935 

1914 

1981 

2030 

12566 

No.  of  events  defined 
by  the  analyst 

2084 

1664 

1495 

1542 

1687 

1754 

10226 

No.  of  events 
accepted  without 
modifications 

1549 

1336 

1267 

1287 

1378 

1354 

_ 

8171 

Table  3.6.1.  IMS  phase  detections  and  event  summary. 


3.7  ESAL  time-lag  statistics 

During  two  years  of  continuous  operation,  the  operational  performance  of  the  IMS  Expert 
System  for  Association  and  Location  (ESAL)  has  in  general  been  very  good.  When  data 
are  available  from  all  arrays,  and  all  EP-SigPro  processes  are  working  normally,  ESAL 
processes  phase  detections  at  close  to  real-time  speed. 

Prom  time  to  time,  however,  some  process  that  provides  ESAL  with  input  data  will  fail  or 
slow  down  and  ESAL  will  fall  behind  real-time,  waiting  for  the  detections  to  become 
available.  The  Chinese  test  explosion  on  5  October  1993,  for  example,  was  not  defined 
until  8  hours  after  the  event  occurred  because  the  processing  of  ARCESS  detections  was 
behind. 

To  monitor  to  what  degree  ESAL  keeps  up  with  real-time,  recording  of  time-lag  was 
started. 

Every  ten  minutes,  the  UNIX  system  time  (based  on  a  GPS  satellite  clock)  is  written  to  a 
file  together  with  the  detection  time  last  processed  by  ESAL  (found  in  the  timestamp  data¬ 
base  table).  Once  a  week,  the  results  are  plotted  and  for  each  peak  in  the  plot,  an  explana¬ 
tion  is  sought. 

The  first  three  of  these  weekly  plots  are  displayed  in  Figs.  3.7.1  -  3.7.3 

The  sawtooth  appearance  of  the  plots  is  a  result  of  die  way  phase  detections  are  processed 
in  the  IMS:  Each  array  has  a  separate  EP-SigPro  process  which  extracts  detection  parame¬ 
ters  and  writes  than  to  a  SIGPRO  database  common  to  all  arrays.  Once  an  hour,  the 
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database  used  by  BSAL.  ESAL  wakes  up  every  30  minutes,  reads  the  new  phase  detec¬ 
tions,  defines  events  from  them  and  goes  back  to  sleep  for  another  30  minutes.  There  is 
therefore  an  inherent  delay  of  up  to  90  minutes  in  the  processing  scheme.  In  addition, 
some  time  is  spent  by  EP-SigPro  processing  each  detection  before  it  is  written  to  the  data¬ 
base  and  lastly,  ESAL  processes  all  die  detections  in  the  batch  before  it  updates  the  times¬ 
tamp  database  table.  As  can  be  seen  in  the  plots  die  combined  effect  is  that  ESAL  during 
normal  operation  stays  between  approximately  30  minutes  and  2  hours  behind  real-time. 

The  plots  also  show  examples  of  several  different  reasons  for  time-lag  peaks:  Computer 
hardware  errors,  missing  data,  too  many  detections  and  software  stops.  In  the  following,  a 
detailed  explanation  of  the  reason  for  each  peak  is  provided 

Oct  12-13:  The  number  of  ARCESS  phase  detections  rose  from  20-50  per  hour  to  150-200 
for  a  4  hour  period.  This  was  probably  caused  by  a  temperature-fall  and  ice-cracking.  The 
extra  processing  load  was  too  much  for  the  computer  running  the  Ep-SigPro  process  and 
ESAL  had  to  wait  for  the  detections  to  become  available. 

Oct  14-16:  Too  many  detections  at  the  Apatity  array.  The  extra,  top  peak  on  Oct  15  came 
from  a  hardware-error  on  an  IMS  computer  (a  broken  CPU  on  njaard).  All  IMS  processes 
were  shut  down  while  the  machine  was  being  repaired. 

Oct  17:  No  data  available  from  GERES S.  The  transmission  line  from  Germany  was  down. 
Since  GERESS  data  are  also  recorded  in  Bochum,  the  missing  data  could  be  retrieved 
when  the  line  came  back  up.  For  arrays  without  local  recording,  data  are  lost  totally  when 
the  line  goes  down.  This  is  apparent  in  the  ESAL  time-lag  plots:  When  the  data  transmis¬ 
sion  comes  back  (Octl7-Octl8)  the  delayed  data  has  to  be  processed  and  there  is  a  slowly 
descending,  saw-toothed  line.  If  ESAL  had  been  waiting  for  data  from  an  array  with  no 
local  recording,  there  would  be  no  data  to  process  and  ESAL  could  have  proceeded  with 
the  other  already-processed  arrays.  (This  is  seen  as  a  straight,  sharp  drop  as  on  Oct  28 
when  ARCESS  came  bade  after  an  outage.) 

Oct  18-19:  No  data  from  the  Apatity  array  because  of  problems  related  to  positioning  the 
satellite  antenna  in  Apatity.  To  bring  the  ES  AL-processing  forward,  it  was  dedded  on  Oct 
19  to  leave  Apatity  temporarily  out  of  the  IMS-processing.  This  accounts  fra  the  sharp 
drop  on  Oct  19. 

The  second,  smaller  peak  late  on  Oct  19  is  due  to  a  stop  in  the  ARCESS  EP-SigPro  pro¬ 
cessing.  The  reason  for  the  stop  is  not  clear,  but  may  have  been  caused  by  bad  data. 

Oct  20:  The  GetArrivals  program  failed  on  two  consecutive  invocations,  no  detections 
were  transferred  from  the  SIGPRO-  to  the  IMS -database  for  three  hours.  The  reason  is 
unknown  but  probably  database-related. 

Oct  21-23:  The  satellite  antenna  problem  in  Apatity  had  been  fixed.  It  was  decided  to 
include  the  Apatity  array  in  the  IMS  processing  again.  ESAL  moved  slowly  forward  while 
die  missing  data  were  being  retrieved  from  the  local  disk-loop  in  Apatity  and  processed. 

Oct  24:  A  full  day  of  continuous,  error-free  processing. 
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Oct  24:  A  full  day  of  continuous,  error-free  processing. 

Oct  25:  ESAL  stopped.  Aborted  in  a  time-interval  where  no  events  woe  defined. 

Oct  26:  Delayed  GERESS  data.  Data-transmission  from  Germany  down. 

Oct  27-28:  ARCESS  data  acquisition  computer  hardware  error.  (Faulty  disk-controller  on 
rein).  The  data  that  were  not  recorded  during  repair  were  lost.  No  processing  required. 

Oct  28-30:  Large  number  of  Apatity  detections  due  to  local  disturbances  and  radio  link 
problems. 

Oct  31-Nov  1 : Close  to  ideal  processing.  No  major  delays. 

This  exercise  of  closely  monitoring  the  ESAL  time-lags  provides  considerable  insight  into 
the  complexity  of  the  task  of  keeping  a  system  like  IMS  running  at  close  to  real-time.  Con¬ 
sidering  the  experimental  and  research-oriented  role  of  die  IMS  operation,  the  priority  is, 
for  die  time  being,  on  collecting  data  from  all  arrays  before  doing  the  network  association 
and  event  location.  This  strategy  may,  however,  change  with  changing  operational  require¬ 
ments. 

U.  Baadshaug 

3.8  GBF  operation 

Events  automatically  located  by  GBF 

The  automatic  GBF  processing  was  temporarily  discontinued  for  a  period  of  47  days  dur¬ 
ing  the  reporting  period,  but  was  otherwise  run  on  a  continuous  basis. 

During  days  091  through  180,  and  222  through  267, 1993, 8445  local  and  regional  events 
were  located  by  GBF.  This  gives  an  average  of  62.0  events  per  processed  day  (136  days 
processed).  69%  of  these  events  are  within  300  km  of  the  nearest  station,  and  84%  of  these 
events  are  within  1000  km  of  the  nearest  station. 

72.2%  of  these  events  were  defined  by  2  regional  phases.  Teleseismic  phases  are  currently 
not  used  by  GBF.  88.8%  of  all  events  had  3  defining  phases  or  less. 


T.  Kvserna 
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4  Improvements  and  Modifications 

4.1  NORSAR 

NORSAR  data  acquisition 

No  modification  has  been  made  to  the  NORSAR  data  acquisition  system. 

The  data  are  recorded  on  a  30-hour  circular  disk  buffer  on  the  IBM  system,  and  archived 
onto  1/2  inch  magnetic  tapes'.  In  addition  to  this,  the  data  are  now  regularly  transmitted  to 
a  SUN  system  for  recording  on  a  48-hour  circular  disk  buffer. 

The  data  from  the  7  subarrays  t  <•  transmitted  over  2400  bps  leased  lines  between  the 
Short  and  Long  Period  Electronic  Modules  (SLEMs)  in  the  subarray  vaults  and  a  Mod- 
comp  computer  at  NDPC.  Severe  problems  have  occurred  with  the  communications,  and 
have  led  to  problems  with  data  quality  and  timing.  An  independent  process  has  been 
developed  which  samples  a  time  pulse  to  automatically  control  the  timing  of  the  old  NOR¬ 
SAR  system. 

The  transmission  between  IBM  and  SUN  has  also  experienced  problems  due  to  old  trans¬ 
mission  devices.  Thus  all  new  processes  developed  for  NORSAR  data  may  be  delayed  due 
to  malfunctions  of  the  IBM-to-SUN  communication  processor. 

NORSAR  detection  processing 

The  NORSAR  detection  processor  has  been  running  satisfactorily  on  the  IBM  4381  com¬ 
puter  during  this  reporting  period. 

Detection  statistics  are  given  in  section  2.  La  addition  to  the  detection  processing  done  on 
IBM,  the  DP  program  is  doing  regular  detection  processing  on  a  SUN  system,  using  the 
unix-based  circular  disk  buffer  (see  below).  A  detection  SNR  threshold  of  20.0  triggers 
automatic  saving  of  waveforms  into  CSS  3.0  data  files. 

NORSAR  event  processing 

There  have  been  no  changes  in  the  routine  processing  of  NORSAR  events,  using  the  IBM 
system. 

In  parallel  with  the  IBM  processing,  routine  event  processing  is  also  done  on  a  SUN  com¬ 
puter  using  the  “old”  IBM  time  delay  correction  data  base  that  has  been  converted  to 
SUN/UNIX.  The  automatic  solutions  produced  are  equal  to  or  better  than  the  old  system 
with  a  lower  false  alarm  rate.  Alert  messages  are  sent  to  USGS  for  events  above  magni¬ 
tude  5.5. 

NORSAR  refurbishment 

We  have  purchased  6  AIM24-1  digitizers  and  one  AIM24-3  digitizer  from  Science  Hori¬ 
zons.  These  7  digitizers  together  with  7  GPS  clocks  have  been  installed  in  the  subarray 
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06 C  vault.  We  are  successfully  recording  data  at  NDPC  using  an  experimental  VSAT 
communication  link.  Details  on  the  refurbishment  effort  are  reported  on  separately. 

Both  short  period  and  long  period  instruments  of  the  current  system  are  functioning  well, 
although  tlx  long  period  instruments  require  a  high  amount  of  manual  calibration.  The 
moat  critical  part  of  the  old  NORSAR  system  is  die  data  acquisition  system,  in  particular 
the  disks,  the  Modcomp  communication  processor  and  the  telecommunication  lines.  The 
refurbishment  project  was  delayed  due  to  lack  of  suitable  digitizers,  but  it  now  seems  to  be 
possible  to  operate  the  AIM24  digitizers  at  the  remote  sites  of  the  NORSAR  array.  After  a 
long-term  field  test,  we  may  be  ready  to  start  installation  of  a  new  data  acquisition  system 
starting  summer  1994. 

42  Regional  Arrays 

Detection  processing 

The  routine  detection  processing  of  the  arrays  is  running  satisfactorily  on  each  of  the 
array’s  SUN-3/280  or  Spares tation  1  acquisition  systems.  The  same  program  is  used  for 
NORSAR,  NORESS,  ARCESS,  FINES  A,  GERESS,  the  Apatity  and  Spitsbergen  arrays, 
but  with  different  “recipes”.  The  beam  table  for  NORESS  and  ARCESS  is  found  in  NOR¬ 
SAR  Sci.  Rep.  No.  1-89/90.  The  beam  table  for  FINESA  and  GERESS  is  found  in  NOR¬ 
SAR  Sci.  Rep.  No.  1-90/91.  The  beam  table  for  the  Apatity  array  is  found  in  NORSAR 
Sci.  Rep.  No.  1-92/93,  and  that  for  the  Spitsbergen  array  is  found  in  NORSAR  Sci.Rep. 
No.  2-92/93. 


Detection  statistics  are  summarized  in  section  3. 

Signal  processing.  Phase  estimation 

This  process  performs  f-k  and  polarization  analysis  for  each  detection  to  determine  phase 
velocity,  azimuth  and  type  of  phase,  and  die  results  are  put  into  the  ORACLE  detection 
and  arrival  tables  for  use  by  the  IMS. 

Event  Processing.  Plot  and  epicenter  determination 

A  description  of  single-array  event  processing  is  found  in  NORSAR  Sci.  Rep.  No.  2-88/ 
89.  and  NORSAR  Sci.  Rep.  No.  2-89/90. 


J.  Fyen 
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5  Maintenance  Activities 

5.1  Activities  in  the  field  and  at  the  Maintenance  Center 

This  section  summarizes  the  activities  at  the  Maintenance  Center  (NMC)  Hamar,  and 
NDPC  activities  related  to  monitoring  and  control  of  NORSAR,  including  monitoring  of 
NORESS,  ARCESS,  FINES  A,  GERES  S  and  the  Apatity  and  Spitsbergen  arrays. 

Activities  involve  preventive  and  corrective  maintenance,  planning  and  activities  related 
to  NORSAR  refurbishment  (NMC/06C).  Preparation  of  a  miniarray  in  Amderma,  Russia, 
involved  assembling  a  6-channel  array  (with  NMC-produced  amplifiers).  In  August  1993 
P.W.  Larsen  (NMC)  joined  a  delegation  to  Russia,  where  they  visited  Dubna  near  Moscow 
and  Peleduy  in  Siberia  in  connection  with  array  siting  surveys. 

The  Spitsbergen  installation  was  completed  in  September  1993. 

NORSAR 

Visits  to  subarrays  in  connection  with: 

•  Replacement  of  RA-S  card 

•  Replacement  of  MP  motors 

•  Replacement  of  FP  motors 

•  Adjustment  of  LP  seismometers,  VE/NS/EW 

•  Adjustment  of  gain  SP/LP  channels 

•  Cable  location 

•  Cable  splicing 

•  NTA/Hamar  assistance  (communications  check) 

•  Demounting  six  telemetry  stations 

•  Line/modem  checks 

•  Installation  of  data  acquisition  equipment  in  connections  with  NORSAR  refurbish¬ 
ment 

NMC 

•  NORSAR  refurbishment  preparations 

•  Preparation  of  a  miniarray  for  Amderma  in  Russia,  involving  assembly  of  a  six- 
channel  array  (dispatched  to  Russia  in  June  1993) 

•  NGI  (Norwegian  Geotechncal  Institute)  survey  (PWL)  at  Kjeller  in  connection  with 
preparation  of  a  3-component  noise  study 


NORESS 

•  Adjustment  of  offset  site  D8 

•  Fiber  optical  link  repair 
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Spitsbergen 

•  Installation  completion 


Subarray/ 

area 

Task 

Date 

NORSAR 

No  visits  to  subarray 

April 

NORESS 

Adjustment  of  DC  offset  site  D8 

20  Apr 

NMC 

Preparations  for  the  NORSAR  refurbishment,  involving  test¬ 
ing  of  different  seismometers  and  digitizers. 

April 

NPDC 

Daily  checks  of  the  following  arrays  have  been  carried  out, 
i.e.,  NORSAR,  NORESS.  ARCESS,  FINES  A,  GERES,  Apat¬ 
ity  and  Spitsbergen. 

SP/LP  instruments  have  been  calibrated  (NORSAR).  Free 
Period  (FP)  and  Mass  Pos.  (MP)  were  measured.  Those  out¬ 
side  specifications  adjusted  (when  feasible  from  NDPC). 

April 

NORSAR 

03C 

Assisted  the  landowner  in  pointing  out  a  cable  in  the  SP  point 
00  area 

28  May 

NMC 

NORSAR  refurbishment  continued 

May 

NDPC 

Daily  checks  of  the  following  arrays  have  been  carried  out, 
i.e.,  NORSAR,  NORESS,  ARCESS,  FINES  A,  Apatity  and 
Spitsbeigen  (partly). 

SP/LP  instruments  have  been  calibrated  (NORSAR),  excl. 
week  20.  FP/  MP  were  measured.  Those  outside  specifica¬ 
tions  adjusted  (when  feasible  from  NDPC). 

May 

NORSAR 

Replaced  RA-S  amplifier  SP  ch  01 

2  June 

02B 

Replaced  MP-motor  vertical  seismometer 

30  June 

02B 

Replaced  MP/FP  motor  NS  seismometer 

30  June 

ARCESS 

NMC-staff  installed  a  new  air-condition  unit 

4-8  June 

NMC 

The  Amderma  miniarray  prepared  at  the  NMC  dispatched  to 
Russia 

June 

June 
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Subarray/ 

area 

Task 

NDPC 

Daily  checks  of  the  following  arrays  have  been  earned  oat, 

Le..  NORSAR,  NORESS,  ARCESS.  FINES  A,  Apatity  and 
Spitsbergen. 

SP/LP  instruments  were  calibrated  (NORSAR),  excl.  week 

26.  FP/  MP  were  measured.  Those  outside  specifications 
adjusted  (when  feasible  from  NDPC). 

June 

NORSAR 

July 

04 C 

SP03  dead  due  to  a  bad  cable.  Checked  gain  of  the  remaining 
SP  channels.  Adjusted  VE/EW  LP  seismometer.  Not  able  to 
operate  NS  seismometer 

1  July 

01B 

SP02  was  dead  due  to  a  bad  cable.  Adjusted  MP/FP  VE/NS/ 
EW  LP  seismometer 

2  July 

04C 

Spliced  cable  SP03 

5  July 

01B 

Spliced  cable  SP02 

7  July 

02C 

Replaced  RA-5  SP01 

8  July 

02C 

Assisted  NTA/Hamar  in  connection  with  a  communications 
check.  Also  checked  LP-channels 

9  July 

04C 

Replaced  MP/FP  motors  (RCDs),  NS  and  EW  LP  seismo¬ 
meters 

26,27  July 

NMC 

P.W.  Larsen  (NMC)  took  part  in  a  survey  at  Kjeller  28  July  in 
connection  with  preparing  a  3 -component  noise  study  instal¬ 
lation. 

28  July 

NDPC 

Daily  checks  of  the  following  arrays  have  been  carried  out, 
i.e.,  NORSAR,  NORESS,  ARCESS,  FINES  A,  Apatity  and 
Spitsbergen. 

SP/LP  instruments  were  calibrated  (NORSAR),  excl.  weeks 
30, 31 .  FP/  MP  were  measured.  Those  outside  specifications 
adjusted  (when  feasible  from  NDPC). 

July 

NORSAR 

August 

02B(telem.) 

Demounted  the  six  telemetry  stations  and  brought  the  equip¬ 
ment  to  the  NMC 

August 

02B 

Line  check,  Adjusted  gain  SP-channels  1, 2, 4  and  5 

11  Aug 

01A 

Line  check 

31  Aug 

NORESS 

Repaired  fiber  optic  link  D3 

25  Aug 
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Su  barmy/ 


In  cooBctiop  with  NORSAR  refurbishment.  Science  Horizons 
data  acquisition  equipment  was  installed  in  the  06C  CTV. 

August 

P.W.  Larsen  visited  Dubna  near  Moscow  in  connection  with 
plans  for  upgrading  the  existing  Pdeduy  station. 

9-20  Aug 

Daily  checks  of  the  following  arrays  have  been  carried  out, 

Le.,  NORSAR,  NORESS,  ARCESS,  FINES  A,  Apatity  and 
Spitsbergen. 

SP/LP  instruments  were  calibrated  (NORSAR),  excL  week 

32..  FP/  MP  were  measured.  Those  outside  specifications 
adjusted  (when  feasible  from  NDPC),  excl.  week  32. 

August 

September 

Repaired  SP04  cable 

16,17,27. 

29  Sept 

Repaired  SP03  dable 

16,17,27, 

29  Sept 

Replaced  RA-5  SP03 

21  Sept 

Pointed  out  a  cable  for  the  power  company.  Adjusted  ch.  gain 
SP03  and  LP/NS  seismometer.  Adjusted  MP  all  LP  instru¬ 
ments 

30  Sept 

Replaced  Remote  Centering  Device  (RCD),  EW  MP/FP 

14  Sept 

Replaced  damping  coil  EW  seismometer 

15  Sept 

Adjusted  ch.  gain  sp  ch  1, 4  and  LP  vertical  channel 

20  Sept 

Adjusted  MP  and  FP  all  LP  seismometers.  Soldered  bad  con¬ 
tact  on  data  coil  vertical  LP  seismometer 

20  Sept 

Cable  sp03  spliced.  Adjusted  ch.  gain  sp  ch  1 , 3  and  all  LP 
channels.  Adjusted  MP  and  FP  all  seismometers 

10  Sept 

Repaired  fiber  optical  link  A2,  A3,  B2,  B4,  C2,  D6  and  D7 

7-9  Sept 

Installation  completed.  Cables  have  been  laid  down  to  the 
remote  sites  Al,  A2  and  B3.  Only  sites  AO,  A2,  B2,  B4  and 

B5  have  good  instruments.  The  remaining  were  damaged  dur- 
ing  transport 

1-6  Sept 

NORSAR 

01A 


ARCESS 
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Subarray/ 

area 

Task 

D 

NPDC 

Daily  checks  of  die  following  arrays  have  been  carried  out, 

Le..  NORSAR,  NORESS.  ARCESS,  FINBSA,  Apatity  and 
Spitsbergen. 

SP/LP  instruments  were  calibrated  (NORSAR).  excl.  week 

37.  FP/  MP  were  measured.  Those  outside  specifications 
adjusted  (when  feasitde  from  NDPC),  excl.  week  37. 

September 

Table  5.1.  Activities  in  the  field  and  the  NORSAR  Maintenance  Center,  including  NDPC 
activities  related  to  NORSAR,  NORESS,  ARCESS,  FINES  A,  GERESS  and  the 
Apatity  and  Spitsbergen  arrays  1  April  -  30  September  1993. 


5.2  Array  status 

As  of  30  September  1993  the  following  NORSAR  channels  deviated  from  tolerances: 

01A  01  8  Hz  filter 

02  8  Hz  filter 

04  30  dB  attenuation 

01B  Out  of  operation 

02B  07 
08 
09 

02C  Out  of  operation 

03C  04 
08 

04C  08 

06C  Subarray  out  of  operation  from  27  September. 


O.A.  Hansen 
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7  Summary  of  Technical  Reports  /  Papers  Published 

7.1  Continuous  threshold  monitoring  of  the  Lop  Nor,  China,  test  site 

Introduction 

The  continuous  threshold  monitoring  technique  (Ringdsl  and  Kverna,  1989)  represents  a 
new  approach  toward  achieving  reliable  seismic  monitoring  for  the  purpose  of  verifying 
nuclear  test  ban  treaties. 

Traditionally,  seismic  monitoring  has  relied  upon  applying  signal  detectors  to  individual 
stations  within  a  monitoring  network,  associating  detected  phases  and  locating  possible 
events  in  the  region  of  interest.  This  procedure  has  been  accompanied  by  assessments  of 
network  capabilities  for  the  target  region,  usually  by  applying  statistical  models  for  the 
noise  level  distribution,  introducing  station  corrections  for  signal  attenuation  and  devising 
a  combinational  procedure  to  determine  the  detection  threshold  as  a  function  of  the  num¬ 
ber  of  phase  detections  required  for  reliable  location. 

The  statistical  noise  models  used  in  these  capability  assessments  are  not  able  to  accommo¬ 
date  the  effect  of  interfering  signals,  such  as  the  coda  of  huge  earthquakes,  which  may 
cause  the  estimated  thresholds  to  be  quite  unrealistic  at  times.  Furthermore,  only  a  statisti¬ 
cal  capability  assessment  is  achieved,  and  no  indication  is  given  as  to  particular  time  inter¬ 
vals  when  die  possibility  of  undetected  clandestine  explosions  is  particularly  high. 

The  continuous  threshold  monitoring  technique  alleviates  these  problems.  It  makes  it  pos¬ 
sible  to  ascertain,  at  any  point  in  time,  for  a  given  target  region,  the  maximum  magnitude 
of  a  possible  clandestine  explosion  at  a  predefined  level  of  confidence.  This  makes  it  pos¬ 
sible  to  focus  attention  upon  those  specific  time  intervals  when  realistic  evasion  opportu¬ 
nities  exist,  while  retaining  confidence  that  no  treaty  violation  has  occurred  at  other  times. 

The  continuous  threshold  monitoring  u.  .  hnique  has  previously  been  applied  experimen¬ 
tally  in  connection  with  the  Novaya  -iemiyfc  test  site  (Ringdal  and  Kvcraa,  1992;  Kvtcrna, 
1992).  This  test  site  is  within  region;/  fiance  of  die  Fennoscandian  arrays,  and  conse¬ 
quently  an  excellent  monitoring  capability  (mb  -  2.5)  has  been  achieved  for  this  site. 

Application  to  the  Lop  Nor  test  site 

In  order  to  further  demonstrate  how  continuous  threshold  monitoring  could  be  performed 
in  a  practical  operation  situation,  we  have  conducted  an  experiment  during  which  we  have 
applied  continuous  threshold  monitoring  to  the  Chines  test  site  at  Lop  Nor  for  a  five-day 
period.  Our  data  base  has  been  the  regional  array  network  NORESS,  ARCESS  and  GER- 
ESS.  As  illustrated  in  Fig.  7.1.1,  these  three  arrays  are  all  at  teleseismic  distances  from  the 
test  site,  with  excellent  P-phase  detection  capabilities  (see  Fig.  7.1.2).  In  particular,  the 
NORESS  array  has  an  excellent  detection  capability  for  this  test  site. 

The  parameters  used  in  the  threshold  monitoring  experiment  are  given  in  Table  7.1.1.  For 
each  array,  we  stem'  “optimum”  P  beams  towards  the  test  site,  and  calibrate  these  beams 
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lining  actually  observed  signal  attenuation  from  previous  Lop  Nor  explosions.  By  focusing 
in  this  way  on  the  target  region,  we  can  at  any  point  in  time  measure  the  “noise  magni¬ 
tude”  for  a  given  phase  at  a  given  array,  and  combine  these  data  to  obtain  a  network 
threshold  as  explained  in  detail  by  Ringdal  and  Kvcraa  (1989). 

Results 

Pigs.  7. 1.3-7. 1.7  show  the  results  of  the  monitoring  experiment.  Each  of  these  figures  cov¬ 
ers  one  data  day,  starting  1  October  1993.  The  upper  three  traces  of  each  figure  represent 
the  thresholds  (Le„  90%  upper  magnitude  limits)  obtained  from  the  three  individual 
arrays,  whereas  the  bottom  trace  illustrates  the  network  threshold.  Typically,  the  individ¬ 
ual  array  traces  have  a  number  of  significant  peaks  for  each  24-hour  period,  due  to  inter¬ 
fering  events  (local  or  teleseismic).  On  the  network  trace,  the  number  and  sizes  of  these 
peaks  are  greatly  reduced,  because  an  interfering  event  will  usually  not  provide  matching 
signals  at  all  the  stations.  Prom  probabilistic  considerations,  it  can  in  such  cases  be 
inferred  that  the  actual  network  threshold  is  lower  than  these  individual  peaks  might  indi¬ 
cate. 

We  will  not  discuss  in  detail  the  individual  peaks  cm  the  network  trace.  Here,  we  will  just 
note  that  on  2  October  (day  27S)  an  aftershock  sequence  occurred  in  the  S.  Sinkiang  prov¬ 
ince.  Furthermore,  the  last  day,  S  October  1993,  was  the  day  of  an  actual  nuclear  explosion 
(tn(,  =  3.9)  at  Lop  Nor,  and  this  event  naturally  stands  out  cm  the  plot.  The  peak  value  of 
the  network  threshold  plot  does  not  represent  the  actual  magnitude  of  die  event,  but  is 
slightly  lower  (see  discussion  below). 

As  a  general  comment  to  Pigs.  7 . 1 .3-7 . 1 .7,  we  note  that  such  plots  are  a  useful  supplement 
to  the  Intelligent  Monitoring  System  (IMS)  (Bache  et  al,  1993),  and  will  enable  the  analyst 
to  obtain  an  instant  assessment  of  the  actual  threshold  level  of  the  monitoring  network. 
The  peaks  on  the  network  traces  may  be  quickly  correlated  with  the  IMS  detection  bulle¬ 
tin,  in  order  to  decide  whether  they  originate  from  interfering  events  or  from  events  in  the 
target  region. 

Discussion 

In  a  monitoring  situation,  it  will  be  important  to  isolate  and  analyze  more  extensively 
those  time  intervals  which  offer  significant  evasion  opportunities.  Table  7.1.2  gives  a  sta¬ 
tistic  of  the  number  of  occasions  during  which  the  upper  magnitude  limit  exceeded  a  given 
level.  In  theory,  if  this  limit  is,  e.g.,  at  4.0,  it  might  be  possible  that  a  clandestine  mj,  =  4.0 
explosion  had  occurred  without  being  detected.  There  are  many  options  available  to  inves¬ 
tigate  such  a  hypothesis  in  more  detail,  although  we  have  not  attempted  to  do  so  in  this 
study.  The  most  immediate  approach  would  be  to  analyze  high-frequency  signals  for  the 
time  interval  being  considered. 

It  is  significant  that  the  3-array  network  studied  in  fills  paper  can  monitor  the  Lop  Nor  test 
site  down  to  3.5  or  below  more  than  99%  of  the  time  (Fig.  7.1.8).  Further  improve¬ 
ments  would  clearly  be  possible  by  adding  more  stations  to  the  monitoring  network,  espe¬ 
cially  highly  sensitive  stations  at  other  azimuths  than  those  covered  by  the  northern 
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European  network.  This  would  in  particular  contribute  to  lowering  the  peaks  due  to  inter¬ 
fering  events,  whereas  any  event  truly  originating  in  the  target  region  would  of  course  still 
stand  out  clearly  on  the  combined  network  traces. 

As  a  final  comment,  we  will  address  the  apparent  contradiction  that  the  magnitude  at  the 
time  of  the  nuclear  explosion  on  the  threshold  trace  is  slightly  lower  than  the  actual  event 
magnitude.  In  order  to  explain  this,  we  recall  that  the  network  TM  calculation  assumes 
that  a  “hidden”  signal  is  “leas  than  or  equal  to”  the  actually  observed  trace  value  for  each 
station.  While  this  is  a  correction  assertion  for  each  individual  station,  it  can  create  a  bias 
if  used  in  a  network  context,  assuming  that  there  is  a  detectable  signal  present. 

Strictly  speaking,  this  model  should  only  be  used  during  periods  with  non -detectable  sig¬ 
nals,  or  when  detections  occur  from  events  outside  the  target  area.  If  there  is  a  detection 
that  could  possibly  correspond  to  an  event  in  the  target  region,  “equal”  should  be  used 
instead  of  “less  than  or  equal”.  Hence,  the  “worst  case”  (upper  90%  limit)  magnitude  in 
this  case  would  be  the  90  per  cent  quantile  of  the  distribution  for  the  maximum-likelihood 
mt,  estimate. 

If  this  philosophy  is  adopted  in  calculating  the  threshold  traces,  it  will  result  in  a  slightly 
increased  height  of  the  peaks  that  are  consistent  (in  azimuth  and  velocity)  with  events  in 
the  target  region.  The  “background”  threshold  level  will  not  change,  and  the  peaks  that  can 
be  confidently  assigned  to  events  in  other  regions  will  be  reduced  in  the  same  way  as 
before.  The  resulting  threshold  trace  computations  will  be  slightly  more  complex  in  those 
cases  where  peaks  at  individual  stations  occur. 

The  above  considerations  amplify  the  importance  of  using  TM  in  combination  with  a  con¬ 
ventional  detection/location  system.  Used  in  this  way,  a  detectable  event  will  be  processed 
in  the  conventional  way,  whereas  upper  magnitude  limits  of  non-detectable  events  will  be 
provided  by  the  TM  method. 

F.  Ringdal 
T.  Kvseraa 
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Station 

H 

TV. 

Tin* 

App. 

VeL 

Azim. 

Filter 

Config. 

STA_len. 

Tun. 

Tot. 

Sta_calib. 

ARC 

n 

479.7 

13.1 

84.8 

2-4  Hz 

A0.C.D 

2.0 

2.0 

1.908 

GER 

D 

550.5 

16.1 

67.1 

0.8-2. 8Hz 

AO.CJ) 

2.0 

2.0 

2.023 

NRS 

a 

530.5 

19.1 

78.1 

1-3  Hz 

A0.CJ5 

1.0 

2.0 

1.133 

Hr.  time 
App.  vel. 
Azim. 
Filter 
Config. 

STAJen. 
Tim.  tol. 
STA_calib. 


Travel  time  of  phase 

Apparent  velocity  from  broadband  F-K  measurement 
Azimuth  from  broadband  F-K  mesurcment 
Cutoffs  of  bandpass  filter  (3rd  order  Butterworth) 

Array  configuration  used  in  beamforming.  AO.B.C  means 
AOZ,  B-ring  and  C-ring 
STA  length  in  seconds 

Time  tolerance  when  searching  for  maximum  STA 
Calibration  factor  used  when  converting  STA  values 
(in  quantum  units)  to  magnitude 
Magnitude  *  loglO(STA)  +  STA_calib. 


Table  7.1.1,  Parameters  used  in  threshold  monitoring  experiment. 


Day-of-Year 

Total 

274 

275 

276 

277 

278 

mt,  £5.0 

0 

2 

D 

0 

1 

3 

mb  £  4.5 

0 

3 

0 

D 

1 

m 

mb  £4.0 

0 

8 

0 

1 

1 

10 

mb  £  3.75 

1 

13 

0 

2 

3 

19 

Table  7.1.2.  Statistics  of  peaks  in  the  network  threshold  traces. 
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Fig.  7.1.1.  Map  showing  the  location  of  the  Lop  Not  test  site  and  the  three  arrays  used  in 
the  monitoring  experiment. 
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F*g>  7.1.2.  NORESS,  ARCESS  and  GERESS  recordings  of  the  Lop  Nor  explosion  of  5 
Oct  1993  (mb»  5.9). 
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Lop  Nor,  China 


Hout»  (GMT) 
October  1,1993 


Fig.  7.13.  Threshold  monitoring  of  the  Novaya  Zemlya  test  site  for  day  274  (1  October 
1993).  The  top  three  traces  represent  thresholds  (upper  90  per  cent  magnitude  limits) 
obtained  from  each  of  the  three  arrays  (ARCESS,  NORESS,  GERESS),  whereas  die 
bottom  trace  shows  the  combined  network  thresholds. 
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November  1993 


Hours  (GMT) 
October  2, 1993 


Fig.  7.1.4.  Same  as  Fig.  7.1.3,  but  for  day  275  (2  October  1993).  The  large  number  of 
threshold  peaks  are  caused  by  an  earthquake  sequence  in  Sinkiang,  China. 
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Fig.  7.1.6.  Same  as  Fig.  7.1.3,  but  for  da 
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(4  October  1993). 
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Threshold  magnitudes  -  Lop  Nor 


Fig.  7.1.8.  Cumulative  statistics  of  the  network  threshold  magnitudes  for  the  five-day 
period,  1-5  October  1993. 
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7 2  The  Lop  Nor  nuclear  explosion  of  5  October  1993 

Introduction 

As  described  in  Section  7. 1,  we  conducted  a  threshold  monitoring  experiment  for  the  Lop 
Nor  test  site  in  Southern  Sinkiang,  China,  for  several  days  prior  to  the  5  October  nuclear 
explosion.  A  large  earthquake  and  aftershock  sequence  occurred  in  the  Southern  Sinkiang 
province  only  three  days  before  the  explosion,  and  this  gave  an  opportunity  to  make  some 
interesting  comparisons,  elaborated  upon  in  this  study. 

The  Lop  Nor  nuclear  explosion  of  5  October  1993 

The  explosion  took  place  on  5  October  1993,  with  origin  time  02.00  GMT.  Table  7.2.1 
lists  the  basic  parameters  of  the  event  as  provided  by  various  sources.  The  mb  magnitudes 
range  from  5.65  to  5.90.  The  most  accurate  location  is  provided  by  the  PDG  bulletin, 
which  uses  a  world- wide  network  for  location  purposes.  The  solutions  by  the  Intelligent 
Monitoring  System  (IMS)  (Bache  et  al,  1993),  both  automatic  (IMS)  and  after  analyst  pro¬ 
cessing  (ARS),  are  also  listed.  The  NORSAR  automatic  and  reprocessed  solutions  (based 
on  the  new  S UN-based  processing  system)  are  included  in  the  table.  The  S UN-based 
NORSAR  trace  plot  is  shown  in  Fig.  7.2.1. 

Figs.  7.2.2  and  7.2.3  show  plots  of  the  interactive  IMS  processing  results.  The  trace  plots 
of  Fig.  7.2.3  are  based  on  array  beams  for  the  four  arrays  Apatity,  ARCESS,  NORESS, 
GERESS  and  a  single  channel  (AO)  for  Spitsbergen.  The  FIN  ESS  array  was  not  yet  opera¬ 
tional  at  the  time  of  the  explosion. 

Table  7.2.2  summarizes  the  automatic  processing  results  for  the  six  arrays.  The  NORESS 
array  has  the  best  signal-to-noise  ratio  (1376.2)  for  this  event,  and  by  extrapolation  this 
array  would  be  expected  to  have  a  detectable  signal  for  an  event  about  2.5  magnitude  units 
lower.  ARCESS  and  NORSAR  also  show  outstanding  SNR.  The  velocity/azimuth  esti¬ 
mates  are  within  the  expected  uncertainty  for  all  arrays.  Note  that  the  Spitsbeigen  array  is 
only  partially  installed,  and  this  is  reflected  in  its  processing  results. 

Comparison  with  previous  events 

We  now  proceed  to  make  a  brief  comparison  between  the  5  October  1993  explosion,  the 
large  21  May  1992  explosion  at  the  same  site  and  the  two  larges ts  events  in  the  2  October 
S.  Sinkiang  earthquake  sequence. 

Table  7.2.3  summarizes  the  PDE  parameters  for  these  four  events.  The  21  May  1992 
explosion  is  comparable  in  mb  to  the  main  shock  of  2  October  1993,  and  the  5  October 
1993  explosion  is  comparable  to  the  2  October  1993  aftershock.  This  similarity  is  illus¬ 
trated  by  Fig.  7.2.3,  which  shows  the  NORESS  P-wave  recordings  (AOZ  seismometer)  for 
the  four  events,  all  plotted  to  the  same  scale. 

Fig.  7.2.4  show  long-period  recordings,  from  the  NORESS  broad-band  seismometers  for 
the  four  events.  Again,  the  same  scale  is  used  in  all  four  cases.  The  surface  waves  of  the 
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|  main  earthquake  have  been  “clipped”  (for  display  purposes).  This  figure  is  very  illustra¬ 

tive,  and  the  following  observations  may  be  made. 

•  As  expected,  the  main  earthquake  and  die  large  May  1992  explosion  have  vastly  dif¬ 
ferent  size  of  surface  waves,  in  spite  of  their  similar  m*,  value.  Thus,  discrimination 
based  on  M^mj,  is  simple  in  these  cases. 

•  The  October  1993  Lop  Nor  explosion  can  likewise  be  readily  identified  as  an  explo¬ 
sion  on  the  basis  of  M^ni),  at  NORESS,  either  by  measuring  Mj  on  the  “marginal” 
surface  wave  shown  on  Fig.  7.2.4,  or  by  using  “negative  evidence”  in  the  case  of 
NORESS.  In  fact,  the  long-period  “noise  magnitude”  is  well  below  the  expected 
value  far  any  earthquake  of  corresponding  m(,  value. 

•  The  surface  waves  of  the  2  October  aftershock  cannot  be  measured  cm  NORESS 
recordings,  and  this  event  cannot  be  identified  as  an  earthquake  from  NORESS  data 
using  Mgim^.  The  reason  is  the  large  coda  level  even  one  full  hour  after  the  main 
shock. 

It  has  previously  been  found,  on  the  basis  of  the  GSETT-2  experiment  (see  section  7.6) 
that  a  modern  network  dominated  by  high-frequency  arrays  is  very  efficient  in  detecting 
aftershocks  closely  following  huge  earthquakes  (see  also  Ringdal,  1992).  The  reason  is 
that  the  high-frequency  coda  drops  very  raoidly  after  the  initial  P  onset,  and  the  high-fre¬ 
quency  arrays  are  able  to  exploit  this  drop  c  in  the  detection  processing.  On  the  other 
hand,  the  long-period  coda  stays  at  a  high  level  for  many  hours  following  huge  earth¬ 
quakes,  and  no  efficient  methods  have  been  found  so  far  to  suppress  this  coda  sufficiently 
to  extract  very  small  surface  waves.  This  is  only  one  of  many  examples  illustrating  that  the 
progress  in  recent  years  in  seismic  event  detection  has  not  been  matched  by  a  similar 
progress  in  event  identification. 

J.  Fyen 
F.  Ringdal 
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Ref. 

Origin  date 

Lat 

Loo 

mb 

IMS  (automatic) 

01.59.59.7 

41.386 

89.619 

5.65 

ARS 

01.59.54.8 

41.110 

89.371 

5.65 

NORSAR  SUN  (automatic) 

02.00.01.0 

42.449 

89.195 

5.88 

NORSAR  Rerun  SUN 

02.00.01.3 

41.365 

87.339 

5.83 

PDE 

01.59.56.5 

41.647 

88.681 

5.90 

Table  7.2.1,  Location  estimates  by  various  systems  of  the  5  October  1993  Lop  Nor 

nuclear  explosion.  Two  of  the  estimates  were  made  automatically  (indicated  in  the 
table). 


Array 

Onset  time 

Res 

STA/LTA 

Vet 

Res 

Azi 

Res 

NORESS 

278:02.08.44.398 

1.71 

1376.2 

16.8 

2.36 

77.5 

1.4 

ARCESS 

278:02.07.54.320 

2.66 

566.7 

143 

0.66 

82.2 

-14.9 

GERESS 

278:02.09.04350 

2.62 

121.9 

17.5 

2.71 

65.8 

-2.2 

Apatity 

278:02.07.28.750 

2.16 

177.3 

14.8 

1.50 

1033 

0.9 

Spitsbergen 

278:02.08.23.900 

335 

37.2 

-6.17 

mm 

-2.9 

NORSAR 

278:02.08.44.800 

1.48 

408.9 

14.5 

0.05 

77.4 

13 

■fable  7.2J,  Automatic  detection  list  for  the  Lop  Nor  nuclear  explosion  05  October  1993. 
The  columns  show  array  name,  automatic  SigFro  onset  time,  onset  residual  relative 
to  PDE  origin  time,  maximum  signal-to-noise  ratio  (STA/LTA),  apparent  velocity 
(km/sec),  residual  in  km/sec,  back-azimuth  in  degrees,  back-azimuth  residual.  All 
residuals  are  relative  to  predictions  using  IASPEI91  tables  and  PDE  origin  time  and 
location. 
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Event 

Ref. 

Origin  time 

Lat 

Loo 

mb 

Ms 

Lop  Nor  92 

PDE 

21  May  92  04.59.57.5 

41.604 

88.813 

6.5 

5*0 

Lop  Nor  93 

PDE 

05  Oct  93  01.59.56.5 

41.647 

88.681 

5.9 

4.8 

Main  shock 

TOE 

02  Oct  93  08.42.32.8 

38.141 

88.638 

6.2 

6.3 

Aftershock 

PDE 

02  Oct  93  09.43.19.5 

38.127 

88.502 

IS 

5.3 

Table  7.2.3.  PDE  parameters  for  four  events  discussed  in  the  text. 
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NORSAR  S cL  Hep  1-93/94 


November  1993 


Fig.  72.2.  Map  showing  the  IMS  solution  (after  analyst  review)  of  5  Oct  93  explosion. 
The  great  circle  path  for  the  detecting  arrays  (based  on  P  and  PeP  estimated  azi¬ 
muths)  are  also  shown. 
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Fig.  7-2-3.  P-phase  waveforms  of  the  5  array  SP  traces  (single  sensor  for  Spitsbergen,  oth¬ 
erwise  array  beams)  for  the  5  Oct  93  explosion. 


Oct  2  ■  CSoO  1993-275:08. 47. 37. 532 


W1Z/93  10:75:41  HOXSM) 


Fig.  7.2.4.  NORESS  P-waves  (AOZ  seismometer)  few  the  four  events  discussed  in  the  text. 
All  traces  are  in  the  same  scale.  Note  that,  for  display  purposes,  two  of  the  traces 
have  been  “clipped”. 
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7 3  Observations  at  NORSAR  of  the  22  September  1993  NPE/CKE 

explosion  at  the  Nevada  Test  Site 

Introduction 

On  22  September  1993,  the  US  Department  of  Energy  detonated  a  one  kiioton  conven¬ 
tional  explosion  at  the  Nevada  test  site  (NTS),  in  an  experiment  named  the  Non-Prolifera¬ 
tion  Experiment  (NPE)  (Springer,  1993).  The  experiment  has  also  been  referred  to  as  the 
Chemical  Kiioton  Experiment  (CKE).  The  experiment  was  conducted  in  the  context  of 
future  CTBT/NPT  monitoring,  and  it  has  the  potential  of  providing  data  useful  for  dis¬ 
crimination  between  chemical  and  nuclear  explosions.  The  NPE  explosion  was  detonated 
in  the  N-tunnel  in  Rainier  Mesa  very  close  to  the  hypocenters  of  several  previous  nuclear 
explosions,  some  of  which  had  yields  similar  to  that  of  the  NPE  explosion.  Through  com¬ 
parison  with  observations  from  previous  shots  at  NTS,  the  NPE  experiment  should  thus 
offer  an  opportunity  to  determine  differences  in  characteristics  of  chemical  and  nuclear 
explosions. 

The  NPE  explosion  was  recorded  on  many  especially-deployed  instruments  in  the  local 
and  regional  distance  range.  Preliminary  reports  indicate  that  the  shot  was  recorded  with 
amplitudes  that  were  in  general  larger  than  those  that  would  be  expected  from  a  similar- 
size  nuclear  shot. 

The  NPE  explosion  was  automatically  detected  and  located  by  some  of  the  arrays  in  north¬ 
ern  Europe  contributing  data  to  the  NORSAR  Data  Processing  Center  (NDPC).  The  pur¬ 
pose  of  this  short  contribution  is  to  present  relevant  detection  and  event  location  data  for 
these  stations,  and  to  make  comparisons  with  observations  for  a  nuclear  test  conducted  at 
NTS  in  1992. 

Data  analysis 

The  nuclear  explosion  that  we  will  compare  the  NPE/CKE  with  is  the  explosion  con¬ 
ducted  at  NTS  on  18  September  1992  and  referred  to  as  “Hunters  Trophy”.  According  to 
the  EDR  of  the  USGS,  this  explosion  was  conducted  at  17:00:00.008  GMT,  at 
37°12 ’24.93 ”N,  116°12’35.94”W,  surface  elevation  2239  m  and  depth  of  burial  385  m. 
Magnitude  is  given  as  mj,  4.4.  The  NPE/CKE  explosion  was  conducted  within  1  km  of 
“Hunters  Trophy”  and  at  about  the  same  depth  of  burial. 

At  the  time  of  the  “Hunters  Trophy”  explosion,  NDPC  received  and  processed  automati¬ 
cally  data  from  the  NORSAR  large-aperture  array  as  well  as  the  high-frequency  arrays 
NORESS,  ARCESS,  FINESA  and  GERESS.  At  the  time  of  the  NPE/CKE  event,  data 
from  two  additional  high-frequency  arrays  were  available,  namely,  the  arrays  at  Apatity 
and  Spitsbergen.  The  FINESA  array  was  not  operational  at  the  time  of  the  NPE/CKE 
event,  due  to  work  related  to  a  refurbishment  of  this  array. 

At  NDPC  all  data  are  subject  to  standard  detection  processing  including  beamforming,  fil¬ 
tering  and  estimation  of  STA/LTA  ratios  for  signal  detection  (EP_SigPro  processing). 
Tables  7.3.1  and  7.3.2  show  automatic  detection  parameters  for  the  “Hunters  Trophy”  and 
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NPE/CKE  events,  respectively.  It  is  seen  that  the  events  are  detected  at  NORSAR,  NOR- 
ESS  and  GERESS  (and  “Hunters  Trophy”  at  FINES  A)  at  distances  ranging  from  73  to  83 
degrees.  The  high-frequency  small- aperture  NORESS  array  has  not  only  the  best  SNR, 
but  broadband  frequency  wave  number  analysis  is  capable  of  estimating  back  azimuth  and 
apparent  velocity  about  as  well  as  the  larger  NORSAR  array.  It  should  be  noted,  however, 
that  the  automatic  processing  of  the  NORSAR  array  is  using  the  full  array.  With  site-spe¬ 
cific  processing  using  the  “best”  subarrays,  NORSAR  beam  SNR  can  be  improved.  From 
Figs.  7. 3. 1  and  7.3.2  we  see  that  NORSAR  subarray  NC6  colocated  with  NORESS  has  a 
larger  SNR  than  the  other  subarrays.  However,  for  event  location  the  full  NORSAR  array 
will  be  superior  to  that  of  individual  subarrays.  ARCESS  detected  neither  of  the  two 
events.  The  NPWCKE  event  was  not  detected  on  the  Apatity  nor  the  Spitsbergen  array. 

The  data  from  the  high-frequency  arrays  are  processed  using  the  Intelligent  Monitoring 
System  (IMS).  Arrivals  on  the  various  arrays  are  automatically  associated  to  form  events 
at  both  local,  regional  and  teleseismic  distances.  IMS  results,  as  reviewed  by  the  analyst, 
are  given  in  Figs.  7.3.3  and  7.3.4.  It  can  be  seen  that  the  event  solution  for  “Hunters  Tro¬ 
phy”  is  better  than  that  for  the  NPE/CKE  event.  The  main  reason  appears  to  be  the  fact 
that  FINESA  was  not  operational  at  the  time  of  the  NPE/CKE  event.  The  mj,  magnitudes 
as  determined  by  IMS  are  4. 14  and  4. 10,  for  the  “Hunters  Trophy”  and  NPE/CKE  events, 
respectively. 

Concluding  remarks 

The  NORSAR,  NORESS  and  GERESS  e  Ttys  at  epicentral  distances  ranging  from  73  to 
83  degrees  all  detected  both  the  “Hunters  llrophy”  explosion  and  the  NPE/CKE  event 
automatically.  The  NPE/CKE  event  was  recorded  with  amplitudes  that  were  only  slightly 
smaller  (generally  by  0.1  tty,  units)  than  those  of  “Hunters  Trophy”  (see  e.g.  Figs.  7.3.1 
and  7.3.2,  where  the  scaling  factors  to  the  left  of  the  traces  can  be  directly  compared  for 
the  two  events.  Amplitudes  at  NORESS  can  be  compared  through  the  scaling  factors  in 
Fig.  7.3.5.)  With  an  m^  of  4.4  as  given  by  USGS  for  “Hunters  Trophy”,  it  is  clear  that  also 
NORSAR ’s  recordings  of  the  NPE/CKE  event  confirm  that  this  one  kiloton  conventional 
explosion  produced  amplitudes  somewhat  larger  than  those  expected  from  a  one  kiloton 
fully  contained  nuclear  explosion. 

As  expected,  the  signals  recorded  at  NDPC  for  these  two  events  had  SNRs  that  were  too 
low  to  permit  meaningful  attempts  at  discriminating  between  their  nuclear  and  chemical 
origins  through  analysis  of  signal  characteristics.  This  is  illustrated  in  Fig.  7.3.5,  which 
shows  the  optimum  beams  at  NORESS  for  these  two  events. 

J.  Fyen 
S.  Mykkdtveit 
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Array 

Onset  daw 

Res 

STA/LTA 

Vd 

Res 

Azi 

Res 

NORSAR 

262:17.11.31.631 

-0.83 

4.8 

17.7 

-1.04 

320.9 

2.4 

NORESS 

262:17.11.33.970 

-0.47 

7.9 

17.5 

-1.36 

323.3 

Si 

ARCESS 

No  detection 

FINESA 

262:17.11.53.230 

-1.57 

4.5 

28.8 

9.06 

316.8 

-13.1 

GERESS 

262:17.12.29.348 

-0.21 

5.8 

18.9 

-2.74 

21.8 

59.9 

Apatity 

Not  installed 

Spitsbergen 

Not  installed 

Table  133.  Automatic  detection  information  for  the  “Hunters  Trophy”  event  of  18  Sep¬ 
tember  1992.  The  columns  show  array,  automatic  EP_SigPro  onset  time,  onset  time 
residual,  detection  STA/LTA,  apparent  velocity  in  km/s,  residual  in  km/s,  back  azi¬ 
muth  in  degrees  and  back  azimuth  residual.  All  residuals  are  relative  to  the 
IASPEI91  travel  time  tables  and  USGS  event  parameters. 


Array 

Onset  time 

Res 

STA/LTA 

Vd 

Res 

Azi 

Res 

NORSAR 

265:07.12.32.608 

0.09 

3.5 

17.6 

-1.21 

322.9 

4.6 

NORESS 

265:07.12.34.212 

-.028 

6.7 

15.9 

-2.98 

323.4 

4.8 

ARCESS 

No  detection 

FINESA 

Not  operational 

GERESS 

265:07.13.29.073 

-0.32 

5.2 

13.0 

-8.67 

342.5 

20.7 

Apatity 

No  detection 

Spitsbergen 

No  detection 

Table  132.  Automatic  detection  information  for  the  NPE/CKE  event  of  22  September 
1993.  The  columns  show  array,  automatic  EP_SigPro  onset  time,  onset  time  resid¬ 
ual,  detection  STA/LTA,  apparent  velocity  in  km/s,  residual  in  km/s,  back  azimuth  in 
degrees  and  back  azimuth  residual.  All  residuals  are  relative  to  the  LASPEI91  travel 
time  tables.  The  origin  time  is  assumed  to  be  07:01:00.000,  and  location  is  assumed 
to  be  the  same  as  that  for  “Hunters  TYophy”. 
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C027  19^2-265:17.11.31.631  Snr  1  24  1.49 

Beom  Vel  17.74  Azi  320.90  P  Amp  56.14 

Lot  41.99  Lon  -116.66  Dep  0.0  Otime  1992-262:17.00.25.564 
AZj  320.90  Dst  7658.8  km  68.93  deg  mop  4,05  2.10  3.62 
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Fig.  73.1.  NORSAR  array  automatic  event  plot  for  the  “Hunters  TYophy”  event.  The 
traces  show  from  top  to  bottom  the  unfiltered  (full)  array  beam,  the  filtered  (full) 
array  beam,  followed  by  seven  filtered  subarray  beams  (subarrays  NB2,  NA0,  NB0, 
NC2,  NC3,  NC4  and  NC6).  The  passband  used  is  1.2  to  3.2  Hz.  Note  that  subarray 
NA0  was  out  of  operation  at  the  time  of  this  event.  The  vertical  line  marks  the  auto¬ 
matic  detection  time,  with  reference  to  subarray  NB2.  The  automatic  event  location 
given  is  within  about  4  degrees  of  the  true  location. 
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Fig.  73.2.  NORSAR  array  automatic  event  plot  for  the  NPE/CKE  event  See  caption  for 
Fig.  7.3. 1  for  details  of  figure  content  and  passband.  Note  that  subarrays  NBO  and 
NC2  were  out  of  operation  at  the  time  of  this  event.  The  automatic  event  location  is 
within  about  6  degrees  of  the  true  location. 


Fig.  733.  IMS  results  for  the  “Hunters  Trophy”  event,  as  reviewed  by  an  analyst  Detec¬ 
tions  with  associated  slowness  values  at  NORESS,  FINESA  and  GERESS  are  used 
to  form  the  event  solution. 
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^*8*  7.3.4.  IMS  result  for  the  NPE/CKE  event,  as  reviewed  by  an  analyst.  Only  informs* 
tion  from  NORESS  and  GERESS  was  available  to  form  the  event  solution  shown. 
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Fig.  733.  Optimum  NORESS  array  beams  (with  respect  to  steering  parameters)  filtered 
in  toe  passband  1.2  to  3.2  Hz  for  the  “Hunters  TYophy”  event  (top)  and  file  NPE/ 
CKE  event  (bottom).  The  two  traces  are  aligned  so  that  toe  signal  onset  is  5  s  after 
the  start  time  of  toe  trace.  The  numbers  to  toe  left  of  toe  traces  are  amplitude  scaling 
factors. 
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7.4  A  generic  algorithm  for  accurate  determination  of  P-phase  arrival 
times 

Introduction 

A  precise  estimate  of  the  onset  time  of  seismic  phases  is  needed  to  obtain  an  accurate 
event  location.  To  obtain  very  precise  onset  times  for  all  types  of  seismic  signals,  seismo- 
logical  observatories  around  the  world  mostly  rely  on  the  picks  provided  by  their  human 
analysts.  However,  the  increase  in  the  number  of  seismic  stations  worldwide  has  not  been 
followed  up  by  a  similar  increase  in  the  number  of  analysts.  The  availability  and  opera¬ 
tional  use  of  reliable,  automatic  procedures  therefore  become  more  and  more  important. 

In  the  automatic  detection  and  signal  processing  module  (SigPro)  used  for  processing  the 
regional  array  data  at  NORSAR,  a  two-step  onset  time  algorithm  is  in  use.  This  procedure 
consists  of  first  applying  a  series  of  short-term  to  long-term  average  (STA/LTA)  detectors 
in  parallel  to  a  set  of  filtered  beams.  When  one  or  more  of  the  STA/LTA  detectors  exceed  a 
predefined  threshold,  a  phase  detection  is  declared  and  a  detection  time  is  found.  Subse¬ 
quently,  a  time  domain  phase  timing  algorithm  is  applied  to  the  filtered  beam  with  the 
highest  SNR,  using  the  detection  time  as  the  starting  value.  A  detailed  description  of  this 
algorithm  is  found  in  Mykkeltveit  and  B ungum  (1984). 

These  SigPro  estimates  of  the  onset  times  are  subsequently  used  by  the  automatic  phase 
association  and  event  location  procedure  (ESAL)  of  the  Intelligent  Monitoring  System 
(IMS)  (Bache  et  al,  1993)  to  produce  a  fully  automatic  event  bulletin.  The  IMS  currently 
provides  for  joint  processing  of  data  from  six  arrays  located  in  northern  and  central 
Europe,  see  Fig.  7.4.1.  The  events  in  the  automatic  bulletin  are  finally  reviewed  and  cor¬ 
rected  by  the  analyst  using  the  Analyst  Review  Station  (ARS)  of  the  IMS.  Through  the 
analyst  review  we  have  experienced  that  the  phase  onset  times  often  have  to  be  signifi¬ 
cantly  adjusted.  In  order  to  improve  the  precision  of  the  automatic  event  locations  pro¬ 
vided  by  the  IMS  and  in  order  to  reduce  the  analyst’s  workload,  there  is  therefore  a  strong 
need  to  improve  the  precision  of  the  automatic  onset  time  estimates. 

Autoregressive  modelling  has  been  shown  to  provide  a  useful  tool  in  characterizing  seis¬ 
mic  noise  and  signals.  Tjpstheim  (1975a, b)  applied  such  modelling  to  the  seismic  discrim¬ 
ination  problem.  Takanami  (1991)  used  autoregressive  models  for  onset  time  estimation 
for  microearthquake  networks.  Pisarenko  et  al  (1987)  developed  a  general  autoregressive 
onset  time  estimator,  which  was  further  elaborated  by  Kushnir  et  al  (1990).  In  this  study 
we  will  investigate  the  use  and  performance  of  this  onset  time  estimation  method  when 
applied  in  an  automatic  mode  under  various  types  of  conditions. 

In  this  paper,  we  develop  a  generic  procedure  to  reestimate  the  onsets  of  all  types  of  first- 
arriving  P-phases  using  the  SigPro  onset  estimates  as  a  starting  point.  By  applying  the 
autoregressive  likelihood  technique,  we  have  obtained  automatic  onset  times  of  a  quality 
such  that  70%  of  the  automatic  picks  are  within  0.1  s  of  the  best  manual  pick.  For  die  Sig¬ 
Pro  onset  time  procedure  currently  used  at  NORSAR,  the  corresponding  number  is  28%. 
We  conclude  that  automatic  reestimation  of  first-arriving  P-onsets  using  the  autoregressive 
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likelihood  technique  has  the  potential  of  significantly  reducing  the  retiming  efforts  of  the 
analyst. 

Autoregressive  likelihood  estimation  of  onset  time 

Following  Pisarenko  et  al  (1987)  and  Kushnir  et  al  (1990),  the  autoregressive  likelihood 
algorithm  for  onset  time  estimation  is  based  on  regarding  the  signal  onset  as  the  time  when 
the  statistical  features  of  the  observed  time  series  are  abruptly  changed.  For  each  argument 
t  within  a  predefined  search  interval  (/j  ,t2)  of  length  N,  autoregressive  models  of  the 
observations  within  the  intervals  (rlt t)  and  (x,r2)  are  calculated  by  a  Levinson-Durbin 
procedure.  From  the  variances  oj  and  o2  of  the  autoregressive  model  residuals  of  the  two 
time  intervals,  a  maximum-likelihood  algorithm  is  used  to  calculate  the  likelihood  func¬ 
tion  L(t)  in  accordance  with  the  formula 

H t)  =  [T/nCT1(t)-(N--c)/no2(T)  ]  (1) 

where  the  argument  to  the  maximum  of  L(t  )  defines  the  onset  time  of  the  signal,  see 
Fig.  7.4.2. 

The  algorithm  working  on  single  component  data,  hereafter  denoted  ESTON1,  takes  into 
account  changes  in  both  power  and  frequency  content,  and  it  is  therefore  important  that 
the  broadband  signal  waveforms  are  retained.  This  is  very  different  from  the  onset  time 
estimator  currently  used  in  SigPro,  which  only  exploits  power  differences  within  the  nar¬ 
row  frequency  band  with  the  highest  signal-to-noise  ratio  (SNR).  The  algorithm  working 
on  three  component  data,  hereafter  denoted  ESTON3,  is  in  addition  sensitive  to  changes  in 
the  polarization  characteristics  of  the  three-component  observations.  Following  the  rec¬ 
ommendations  of  Pisarenko  et  al  (1987),  we  have  in  all  our  calculations  used  autoregres 
sive  modelling  of  order  3. 

It  is  noteworthy  that  both  ESTON1  and  EST0N3  require  that  the  search  be  limited  to  a 
relatively  short  time  window.  If  an  initial  event  location  and  origin  time  is  known,  we  can 
determine  the  required  short  time  window  for  the  search.  Alternatively,  the  phase  onsets 
provided  by  SigPro  can  be  used  to  restrict  the  search.  In  any  case,  the  autoregressive  like¬ 
lihood  estimation  of  onset  time  should  be  well-suited  to  a  post-processing  application. 

Generic  application;  retiming  of  first-arriving  P-phases 

We  have  conducted  an  experiment  in  reestimating  the  onset  time  of  all  first-arriving  P- 
phases  defined  in  the  automatic  IMS  bulletin,  using  the  ESTON1  method.  For  a  period  of 
four  days  (September  27  -  30, 1993),  391  first-arriving  P-phases  associated  with  events  in 
the  IMS  bulletin  were  defined.  They  were  distributed  among  all  the  arrays  shown  in 
Fig.  7.4. 1,  and  originated  from  events  at  both  local,  regional  and  teleseismic  distances. 
All  P-phases  were  carefully  retimed  using  an  interactive  signal  processing  package  (EP) 
with  high-resolution  graphics  (Fyen,  1989),  and  about  10%  of  them  were  rejected  due  to 
false  detections  or  erroneous  phase  association,  such  that  350  first-arriving  P-phases 
remained  for  further  analysis  after  this  manual  screening  process.  When  comparing  these 


99 


NOR S AH  ScL  Rep.  1-93/94 


November  1993 


numbers  to  the  general  IMS  performance  (Mykkeltveit  et  al,  1993),  it  appears  that  this 
sample  is  fairly  typical  for  an  operational  situation. 

The  149  P-phases  recorded  during  the  two  first  days  of  the  time  period  were  used  to  tune 
the  implementation  of  ESTON1.  By  comparing  the  differences  between  the  manual  and 
the  SigPro  onset  times,  a  maximum  difference  of  2.8  s  was  observed.  Consequently,  the 
search  interval  to  be  used  by  ESTON1  was  set  to  ±3  s  around  the  SigPro  onset. 

The  different  types  of  P-phases  (Pg,  Pn,  P  and  PKP)  spanned  a  wide  range  of  signal  char¬ 
acteristics  with  respect  to  spectral  content,  complexity,  SNR  and  signature  (impulsive, 
emergent).  From  extensive  testing  of  ESTON1,  we  found  that  in  order  to  successfully  pro¬ 
cess  all  types  of  signals,  we  had  to  identify  the  widest  possible  spectral  band  for  which  the 
signal  had  usable  SNR.  This  was  done  in  the  time  domain  by  estimating  the  maximum 
SNR  within  the  search  interval  in  a  series  of  narrow  passbands.  The  spectral  band  was 
defined  such  that  we  initially  selected  the  narrow  frequency  band  with  the  highest  SNR.  If 
the  neighboring  frequency  bands  had  an  SNR  within  a  factor  of  S  of  the  maximum  and 
also  exceeded  an  SNR  of  4,  the  spectral  band  was  extended  so  as  to  include  this  band  as 
well. 

Our  experiments  also  showed  that  in  order  to  obtain  stable  estimates  of  the  likelihood 
function  L  (r) ,  it  was  important  to  filter  and  decimate  the  data  in  accordance  with  the 
highest  frequency  of  the  signal  spectrum.  For  signals  with  a  high  SNR  (typically  above 
40)  and  a  wide  bandwidth,  no  filtering  or  decimation  was  needed. 

We  found  that  the  onsets  provided  by  ESTON1  were  biased  slightly  late,  and  the  delay 
appeared  to  be  linearly  dependent  on  the  dominant  period  of  the  signal.  By  linear  regres¬ 
sion  of  all  signals  with  SNR  >  6,  the  bias  b  could  be  approximated  by  the  relation 
b  ~  0.38p  where  p  is  the  dominant  period  of  the  signal.  The  flowchart  of  Fig.  7.4.3  out¬ 
lines  the  processing  steps  involved  in  the  reestimation  of  the  arrival  time  of  first-arriving 
P-phases  using  the  ESTON1  method. 

The  201  P-phases  recorded  during  the  last  two  days  of  the  test  period  were  used  to  evalu¬ 
ate  the  new  procedure.  Fig.  7.4.4a  shows  the  difference  between  the  manually  picked 
onsets  and  the  automatic  onsets  from  SigPro  versus  the  highest  SNR  measured  in  any  nar¬ 
row  filter  band.  For  comparison,  Fig.  7.4.4b  shows  the  difference  between  the  manually 
picked  onsets  and  the  automatically  reestimated  onset  times  using  the  ESTON1  method. 
From  comparing  these  two  figures  it  is  apparent  that  the  improvement  when  using 
ESTON1  is  significant  for  all  SNRs. 

To  quantify  the  improvement,  we  have  in  Fig.  7.4.5  plotted  the  percentage  of  the  observa¬ 
tions  within  a  range  of  absolute  time  differences  between  the  automatic  and  the  manual 
picks.  For  SigPro,  50  percent  of  the  automatic  onsets  were  within  0.23  s  of  the  manual 
pick,  whereas  for  ESTON1  the  50  percent  level  (median)  was  as  low  as  0.05  s. 

We  also  divided  the  observations  into  a  teleseismic  and  a  local/regional  data  set.  For  Sig¬ 
Pro,  the  median  time  differences  were  about  equal  for  the  two  data  sets.  For  ESTON1,  the 
median  time  difference  was  slightly  smaller  for  the  local/regional  data  set  than  for  the  tel- 
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eseismic.  This  difference  could  be  due  to  generally  longer  dominant  periods  of  the  tele- 
seismic  P-phases. 

As  expected  and  also  seen  from  Figs.  7.4.4a  and  7.4.4b,  the  precision  of  the  automatic 
onsets  is  best  for  high  SNRs.  By  again  dividing  the  observations  into  two  data  sets,  one 
with  SNR  less  than  or  equal  to  10  and  one  with  SNR  greater  than  10,  we  found  that  SigPro 
had  a  median  difference  of  0.29  s  for  the  low  SNR  data  set  and  0.19  s  for  the  other.  The 
corresponding  numbers  for  ESTON1  were  0. 10  s  and  0.04  s,  respectively. 

The  implications  on  the  analyst's  retiming  efforts  can  be  illustrated  by  the  following 
example:  If  we  assume  that  the  analyst  will  accept  a  maximum  deviation  of  0. 1  s  from  the 
"correct"  manual  pick  without  doing  retiming,  we  can  from  Fig.  7.4.5  see  that  28  percent 
of  the  SigPro  onsets  are  acceptable,  whereas  70  percent  of  the  ESTON1  onsets  are  accept¬ 
able.  Clearly,  automatic  reestimation  of  first-arriving  P-onsets  using  the  algorithm 
described  above  has  the  potential  of  significantly  reducing  the  retiming  efforts  of  the  ana¬ 
lyst. 

Conclusions 

The  results  presented  in  this  study  show  that  very  precise  automatic  estimates  of  phase 
onsets  can  be  obtained  with  the  autoregressive  likelihood  estimation  technique.  Imple¬ 
mentation  of  the  method  requires  that  we  have  available  approximate  estimates  of  the 
phase  arrival,  and  we  have  shown  that  such  approximate  estimates  can  be  obtained  from 
automatic  event  definitions  (phase  association  and  event  location)  by  the  Intelligent  Moni¬ 
toring  System  (IMS).  In  this  way  the  autoregressive  likelihood  estimation  method  can  pro¬ 
vide  phase  onsets  that  match  the  human  precision.  This  has  previously  been  demonstrated 
for  events  from  the  Khibiny  Massif,  by  quantifying  the  uncertainty  of  both  manual  and 
automatic  onset  estimates  of  various  phases  at  the  Apatity  stations  and  at  ARCESS 
(Kvasrna,  1993).  Furthermore,  the  precision  of  the  automatic  phase  picks  shows  very  large 
improvement  in  comparison  to  the  automatic  phase  onsets  from  the  continuous  processing 
providing  input  to  the  IMS. 

We  realize  that  in  order  to  obtain  accurate  event  locations,  precise  onset  time  estimates  are 
necessary,  but  not  sufficient.  If  the  theoretical  travel-time  model  used  in  the  event  location 
deviates  from  the  true  travel-times,  the  accuracy  of  the  event  locations  will  be  reduced. 
Introduction  of  travel-time  corrections  as  well  as  other  aspects  of  accurate  event  location 
are  discussed  by  Kveema  and  Ringdal  (1993). 

During  the  work  with  the  autoregressive  likelihood  estimation  method,  we  have  experi¬ 
enced  that  the  display  of  the  likelihood  functions,  as  illustrated  in  Fig.  7.4.2  can  assist  the 
analyst  in  picking  the  correct  phase  onsets.  In  the  context  of  interactive  analysis  of  seismic 
data,  we  believe  that  the  idea  of  making  such  likelihood  functions  available  to  the  analyst 
should  be  pursued. 

It  is  clear  that  when  estimating  arrival  times  by  the  autoregressive  method,  the  results  for 
specific,  well-calibrated  regions  are  more  precise  than  can  be  obtained  when  the  method  is 
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used  in  a  "generic"  mode.  Efforts  should  be  made  to  extend  the  number  of  well-calibrated 
regions  in  order  to  make  such  optimum  use  of  the  method. 

T.  Kvserna 
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Fig.  7.4.1,  Map  showing  the  locations  of  the  six  regional  arrays  currently  used  by  die 
Intelligent  Monitoring  System  at  the  NORSAR  data  processing  center. 
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7.4.2.  The  top  trace  is  the  likelihood  function  resulting  from  autoregressive  onset  time  estimation  of  the  data  in 
the  bottom  trace.  The  maximum  of  the  likelihood  function  corresponds  to  the  estimated  onset  time. 
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First-arriving  P -phases  from  the  automatic  IMS  bulletin 


New  onset  time 


Fig.  7.4.3.  Flowchart  illustrating  the  processing  steps  involved  in  the  reestimation  of  the 
arrival  time  of  first-arriving  P-phases  using  the  ESTON1  method. 
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-3-2-10123 

SigPro  -  Manual  P-onset  (s) 


-3-2-10123 

ESTON1  -  Manual  P-onset  (s) 


Fig.  7.4.4:  This  figure  show  the  time  difference  between  die  automatic  and  the  manually 
picked  onsets  of  the  201  first-arriving  P-pbases  analyzed  in  this  study  plotted  versus 
file  SNR  of  file  signal. 

a)  shows  the  time  differences  between  the  automatic  onsets  from  SigPro  and  the 
manual  picks.  The  median  absolute  time  difference  is  0.23  s. 

b)  shows  the  time  differences  between  the  reestimated  onsets  from  ESTON1  and  the 
manual  picks.  The  median  absolute  time  difference  is  0. 05  s. 
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Absolute  time  difference  (s) 


Fig.  7.4.5.  These  two  curves  show  the  percentage  of  the  automatic  onsets  within  a  range  of 
absolute  time  differences  from  the  manual  picks.  For  SigPro  (dashed  line).  50  per¬ 
cent  of  the  onsets  are  within  0.23  s  of  the  manual  pick,  whereas  for  ESTON1  (solid 
line)  the  50  percent  level  (median)  is  as  low  as  0.05  s. 
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7.5  Onset  time  estimation  and  location  of  events  in  the  Khibiny 
Massif,  Kola  Peninsula,  using  the  Analyst  Review  Station 

Introduction 

The  technique  of  intelligent  poet-processing  of  seismic  events  proposed  by  Kvcma  and 
Ringdal  (1993)  has  been  shown  to  give  a  substantial  improvement  in  location  accuracy 
when  applied  to  seismic  events  in  the  Khibiny  Massif,  Kola  Peninsula.  In  this  paper,  we 
compare  die  performance  of  the  analyst  using  the  Analyst  Review  Station  (ARS)  to  these 
results.  As  part  of  this  study,  we  estimate  the  uncertainty  in  analyst  time  picks  for  phases 
from  various  regional  arrays  (see  Fig.  7.S.1),  and  we  discuss  the  implication  of  these 
uncertainties  in  terms  of  the  resulting  effect  on  location  accuracy.  An  important  conclu¬ 
sion  inferred  from  this  work  is  that,  in  many  cases,  location  accuracy  does  not  improve 
when  adding  new  phase  readings  and  applying  current  location  programs. 

The  Khibiny  Massif  events 

Six  apatite  mines  are  located  within  an  area  of  about  10  km2  in  the  Khibiny  Massif  on  the 
Kola  Peninsula  of  Russia  (see  Pig.  7.5.2).  A  detailed  description  of  these  mines  and  the 
mining  activity  is  found  in  Mykkeltveit  (1992).  Although  we  have  no  explicit  information 
on  the  exact  sizes  of  these  mines,  interpretation  of  various  maps  suggests  that  the  typical 
size  is  about  1  km2.  The  Kola  Regional  Seismological  Centre  has  since  the  beginning  of 
1991  provided  NORSAR  with  information  on  mining  blasts  in  the  six  Khibiny  Massif 
mines.  Detailed  information  on  the  58  events  used  in  this  study  is  given  in  Kvserna  (1993). 

Kvaema  (1993)  investigated  the  potential  automatic  use  of  an  onset  picker  based  on 
autoregressive  likelihood  estimation.  Both  a  single-component  version  and  a  three-com¬ 
ponent  version  of  this  method  were  tested  cm  data  from  events  located  in  the  Khibiny  Mas¬ 
sif,  recorded  at  the  Apatity  array,  die  Apatity  three-component  station  and  the  ARCESS 
array.  Using  this  method,  he  was  able  to  estimate  onset  times  to  an  accuracy  (standard 
deviation)  of  about  0.05  s  for  P-phases  and  0. 15-0.20  s  for  S-phases.  He  noted  that  these 
accuracies  are  as  good  as  the  best  analyst  picks,  and  considerably  better  than  the  accura¬ 
cies  of  the  current  onset  procedure  used  far  processing  of  regional  array  data  at  NORSAR. 

Estimating  the  precision  of  manual  onset  time  picks 

As  reported  by  Kvasma  (1993),  P  and  S  onsets  at  two  stations  in  Apatity,  APA0  and  APZ9, 
and  the  Pn  onsets  at  ARCESS  were  manually  picked  using  die  interactive  EP  program 
(Fyen,  1989).  Given  the  fact  that  the  characteristics  of  the  Khibiny  Massif  events  were 
known,  the  manual  phase  picking  was  considered  to  be  done  under  "optimum  conditions". 
By  "optimum  conditions"  we  mean  that  the  analyst  utilized  information  on  die  approxi¬ 
mate  phase  arrival  times  and  looked  for  typical  signatures  of  the  different  phases.  He  also 
selected  filters  and  seismometer  components  so  as  to  obtain  the  highest  SNR. 

For  the  purpose  of  the  study  reported  in  this  paper,  all  events  were  reviewed  by  another 
analyst  using  the  Analyst  Review  Station  (ARS)  of  the  IMS.  This  analyst  made  time  picks 
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far  available  p******  for  four  array*  (NORESS,  ARCESS,  FINES  A,  APATITY)  as  well  as 
for  the  APZ9  station.  We  consider  tbe  phase  picks  from  the  ARS  to  be  obtained  under  so- 
called  "operational  conditions"  and  they  may  therefore  be  less  precise  than  those  obtained 
nnHw  "optimum  conditions".  This  is  due  to  the  fact  that  ARS  is  used  as  a  tool  for  routine 
analysis  (i.e..  relatively  short  time  spent  on  each  pick)  of  large  quantities  of  data  and  that 
the  analyst  did  not  have  readily  available  information  on  the  characteristics  of  the  Khibiny 
Massif  events. 

Following  Sereno  (1990),  an  unbiased  estimate  of  the  measurement  variance  is  deter¬ 
mined  from  the  arrival  time  difference  between  two  phase  observations  for  repeated 
events  in  die  same  mine. 

Specifically: 


L  Ei ^,-<4Ui2 

„2  +fI2  _  *-l  »-  1 _ 

1  (Nob,-Nmi„s) 


(1) 


where  ctJ  and  are  the  picking  variance  of  each  phase,  AFo6j  is  the  ith  observation  of 
the  arrival  time  difference  for  the  kth  mine.  <A  To6j>k  is  the  mean  arrival  time  difference 
for  the  kth  mine.  Nobs  is  the  total  number  of  observations  (at  all  mines),  and  Nmines  is  the 
number  of  mines. 

Kvssrna  (1993)  used  formula  (1)  in  various  combinations  to  estimate  standard  deviations 
of  time  picks  for  various  phase  types  and  stations.  He  found  that  the  P-phase  at  APAO 
could  be  picked  with  a  precision  of  o  =  0.04  seconds  when  the  pick  was  made  by  the  ana¬ 
lyst  under  “optimum  conditions”.  In  the  present  study,  we  will  use  these  P-times  (for 
APAO)  as  reference,  and  we  will  assume  that  their  standard  deviation  a2pick  is  0.04.  In  this 
way,  we  can  estimate  °l  rick  directly  from  (1)  for  each  mine,  and  average  these  data  over 
the  six  mines  (using  the  number  of  events  as  a  weighting  factor)  to  obtain  overall  esti¬ 
mates  of  die  uncertainty. 

Results 

The  resulting  estimates  of  the  precision  in  time  picks  by  the  analyst,  using  the  ARS  sta¬ 
tion,  are  presented  in  Table  7 .5.1  and  Figs.  7.5.3-7.S.4. 

Fig.  7.5.3  shows  the  results  for  the  Apatity  array  APAO  and  the  3-component  station 
APZ9.  The  array  has  a  better  precision  for  P  phases  (0.05  versus  0.08),  probably  because 
of  a  far  better  P-wave  SNR  (see  Tkble  7.5.1).  However,  the  3-component  station  has  more 
prerise  S  and  Rg  estimates,  probably  because  of  their  more  impulsive  nature  compared  to 
the  array  recordings.  Note  that  these  secondary  phases  have  a  far  lower  accuracy  in  the 
time  picks  than  the  P -phases.  Also  note  that  for  the  P  and  S  phases  the  ARS  analyst  picks 
are  not  as  precise  as  the  automatic  time  picks  presented  by  Kvsraa  (1993). 
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Pig.  7.5.4  shows  the  results  for  Pn,  Sn  and  Lg  at  the  three  arrays  ARCESS,  FINESA  and 
NORESS.  Pg  for  ARCESS  is  also  shown.  Not  unexpectedly,  Pn  has  the  most  precise 
picks,  followed  by  Sn  and  Lg.  The  ARCESS  Pn  has  by  far  the  best  precision.  This  is  rea¬ 
sonable  in  view  of  the  very  high  SNR  for  these  phases  (Table  7.5.1).  Note  also  that  for 
FINESA  and  NORESS  it  is  possible  to  read  phases  for  only  about  half  of  the  events  or 
fewer. 

The  much  larger  uncertainty  in  FINESA  and  NORESS  P-precision  compared  to  ARCESS 
(about  0.8  versus  0.09  seconds)  is  noteworthy.  Clearly,  the  location  program  ought  to  take 
this  difference  into  account  and  weigh  the  data  accordingly.  At  present  no  routine  mecha¬ 
nism  for  doing  such  weighing  is  applied  in  the  IMS  system,  although  the  option  to  do  so 
exists. 

Location  results 

Figs.  7.5.5-7.5.7  show  plots  of  event  locations  obtained  under  three  different  scenarios. 
All  location  estimates  have  been  made  with  an  assumed  0  km  depth. 

Fig.  7.5.5  shows  the  location  provided  by  the  automatic  IMS  system  with  no  analyst 
review.  All  the  arrival  times  used  here  are  taken  directly  from  the  SigPro  processing,  and 
are  thus  subject  of  significant  uncertainty.  The  median  location  error  is  10.6  km,  which 
must  be  considered  excellent  for  a  fully  automatic  system.  Some  “outliers”  are  due  to 
occasional  erroneous  phase  identification  by  the  automatic  system. 

Fig.  7.5.6  shows  die  results  after  applying  intelligent  post-processing  to  the  ARCESS  and 
Apatity  arrays.  The  median  error  is  now  1.9  km  and  the  worst  case  error  is  5.9  km. 

Fig.  7.5.7  shows  the  results  after  using  the  analyst  (ARS)  reviewed  data  in  the  location 
procedure.  The  median  error  is  3.3  km,  and  the  worst  case  error  is  14.5  km.  These  results 
are  much  better  than  the  automatic  IMS  processing,  but  not  as  good  as  for  the  intelligent 
post-processing. 

As  shown  by  Kvsroa  (1993),  the  ARS  picks  for  the  Apatity  stations  and  ARCESS  are  not 
quite  as  good  as  the  intelligent  post-processing  picks.  A  slight  degradation  in  location 
accuracy  must  therefore  be  expected.  Nevertheless,  we  consider  that  a  main  reason  why 
the  ARS  locations  do  not  match  those  of  the  intelligent  post-processing  is  the  inclusion  of 
NORESS  and  FINESA  readings  in  the  data  base,  without  appropriate  weighting.  We  plan 
to  pursue  this  problem  further  in  the  future. 

F.  Ringdal 
T.  Kvsema 
B.Kr.  H  ok  I  and 
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Table  7.5.1 
are: 

N 

o 

SNR 


APATITY  ARRAY  (APAO) 


N 

a 

SNR 

Pg 

58 

0.051 

57.60 

Lg 

58 

0.389 

8.43 

Rg 

57 

0.494 

8.59 

APATITY  3-COMPONENT  STATION  (APZ9) 


N 

G 

SNR 

Pg 

58 

0.080 

15.20 

Lg 

58 

0.184 

7.79 

Rg 

57 

0.254 

8.89 

.  Basic  data  corresponding  to  Figs.  7.5.3  and  7.5.4.  The  entries  in  the  tables 

:  Number  of  phases  analyzed  (outliers  in  parantheses) 

:  Estimated  standard  deviations  (s)  of  ARS  time  picks 

:  Geometric  average  of  the  linear  signal-to-noise  ratio  (STA/LTA)  of  the  N 
phases.  SNR  of  non-detections  have  been  set  to  3.5  for  P-phases  and  2.5  for 
S-phases. 
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Fig.  7.5.1.  Map  showing  the  locations  of  the  six  regional  arrays  currently  used  by  the 
Intelligent  Monitoring  System  at  the  NORSAR  data  processing  center. 
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Fig.  7.5.2.  In  the  upper  part,  a  large  reference  area  is  shown.  The  location  of  the  ARCESS 
array  is  given  by  a  filled  circle,  and  the  location  of  the  Khibiny  Massif  region  is 
shown.  The  lower  part  shows  a  detailed  picture  of  the  Khibiny  Massif  region.  The 
locations  of  the  six  mining  sites  are  given  by  large  numbers  1-6.  The  Apatity  array 
(APA0)  is  shown  as  a  filled  circle  and  the  three-component  station  (APZ9)  in  the 
town  of  Apatity  is  shown  as  a  large  triangle. 


L 


APZ9  APAO 


NORSAR  Set  Rep.  1-93/94 


November  1993 


ARS  picking  uncertainty  (s) 
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Fig.  7.5.3.  Standard  deviations  of  analyst  time  picks  for  stations  APAO  and  APZ9  using 
the  Analyst  Review  Station  (ARS)  for  the  event  data  base  described  in  the  text. 


NORESS  FINESA  ARCESS 
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Fig.  7.5.4.  Standard  deviations  of  analyst  time  picks  for  the  arrays  ARCESS,  FINESA, 
NORESS  using  the  Analyst  Review  Station  (ARS)  for  the  event  data  base  described 
in  the  text. 
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Expert  system  location  errors 
Median:  10.6  km  Maximum:  85.3  km 
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Fig.  7.5.5.  Location  of  the  six  mining  sites  in  the  Khibiny  Massif  (large  numbers  1-6)  and 
the  locations  of  the  58  reference  events  (small  numbers  1-6)  as  given  by  the  auto¬ 
matic  IMS  processing.  In  the  upper  part,  a  large  reference  area  is  shown,  with  the 
mines  plotted  as  filled  squares.  The  lower  part  shows  a  detailed  picture  for  the  area 
near  the  mines.  The  small  numbers  (1-6)  associated  with  each  event  represent  the 
mine  in  which  die  event  actually  occurred.  The  Apatity  array  is  shown  as  a  filled  cir¬ 
cle  and  the  three-component  station  in  the  town  of  Apatity  is  shown  as  a  filled  trian¬ 
gle. 
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ARCESS  and  Apatity  array  location  errors  (uncalibrated) 
Median:  1.9  km  90%  quantile:  3.6  km  Maximum:  5.9  km 
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Fig.  7.5.6.  Same  as  Fig.  7.5.5,  but  showing  location  results  by  the  automatic  post-process¬ 
ing  method  described  by  Kvserna  and  Ringdal  (1993).  Only  ARCESS  and  the  Apat¬ 
ity  array  are  used. 
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Analyst  reviewed  locations 


Median:  3.3  km  90%  quantile:  9.7  km  Maximum:  14£  km 

31*  32’  33’  94*  36*  36' 


Fig.  7.5.7.  Same  u  Fig.  7.5.5.,  but  showing  location  results  using  ARS  analysis  with  four 
regional  arrays  (ARCESS,  APATITY,  NORESS,  FINES  A). 
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7.6  Generalized  Beamforming  as  a  tool  in  I  DC  processing  of  large 
earthquake  sequences 

Introduction 

Generalized  Beamforming  (Ringdal  and  Kvrema,  1989)  is  a  technique  for  joint  processing 
of  time-aligned  waveforms  from  a  seismic  network.  The  time-alignment  is  made  for  a  grid 
of  beampoints,  and  the  density  and  spatial  coverage  of  the  beam  deployment  can  be  set 
without  any  restrictions. 

The  Generalized  Beamforming  (GBF)  method  has  been  applied  successfully  for  phase 
association  and  event  location,  both  at  regional  distances  (Ringdal  and  Kvserna,  1989; 
Kv«ma,  1990, 1992a)  and  in  a  teleseismic  context  (Taylor  and  Leonard,  1992;  Kvsema, 
1992b).  In  this  paper  we  investigate  the  potential  of  the  GBF  technique  in  achieving  a 
rapid,  preliminary  association  of  phases  for  a  large  aftershock  sequence.  As  is  well  known, 
such  sequences  are  often  problematic  to  process  using  conventional  phase  association 
techniques  since  there  are  so  many  individual  phase  detections  that  the  number  of  possible 
combinations  becomes  very  large. 

The  W.  Caucasus  earthquake  sequence,  April  1991 

On  29  April  1991  a  large  earthquake  (Mj  =7.3)  occurred  in  western  Caucasus,  with  coor¬ 
dinates  42.453N,  43.673E,  h  =  17  km  (NEIC). 

The  earthquake  was  followed  by  a  large  number  of  aftershocks.  According  to  the  cata¬ 
logue  of  Starovoit  et  al  (1992),  114  aftershocks  were  recorded  on  the  day  of  the  main 
shock  (29  April)  and  360  aftershocks  had  been  recorded  by  the  end  of  May. 

The  earthquake  occurred  early  during  the  Group  of  Scientific  Experts  (GSE)  Second  Tech¬ 
nical  Test  (GSETT-2,  main  phase,  see  GSE/CRP/190/Rev.4, 1991),  and  caused  a  consider¬ 
able  load  at  the  National  Data  Centers  (NDCs)  as  well  as  the  four  Experimental 
International  Data  Centers  (EIDCs).  The  day  29  April  was  selected  as  one  of  the  days  for 
which  reprocessing  was  to  be  made  at  EIDCs.  Consequently,  this  day  is  useful  for  study¬ 
ing  the  performance  of  the  experimental  global  system  during  a  day  of  particularly  high 
seismic  activity.  Moreover,  it  provides  an  excellent  opportunity  to  evaluate  the  GBF  tech¬ 
nique  applied  to  a  huge  aftershock  sequence. 

Method 

We  selected  1 1  stations  from  the  total  of  60  participating  in  GSETT-2  for  this  analysis  (see 
Fig.  7.6.1).  These  11  stations  comprised  those  that  had  the  best  detection  performance  for 
the  W.  Caucasus  area.  Table  7.6. 1  lists  the  stations  and  summarizes  the  GBF  parameters 
for  this  experiment.  Note  that  only  one  generalized  beam  was  framed,  and  it  was  steered 
to  42  JN  43 .5E.  The  time  and  azimuth  tolerances  were  set  in  accordance  with  the  GSE 
requirements,  and  adjusted  for  the  beam  focus  area  of  0.3  degrees  radius.  These  tolerances 
were  narrow  enough  to  avoid  many  false  associations,  while  still  allowing  fra  the  typical 
uncertainty  in  detection  times  and  automatic  parameter  estimates.  Detection  threshold  was 
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set  at  3  ™»trhing  phases,  and  GBF  detections  lets  than  IS  seconds  apart  were  grouped 
together. 

Table  7.6.2  shows  the  detection  list  generated  by  the  automatic  GBF  process  for  the  day  in 
question.  For  each  line  our  assessment  of  the  detection  is  given  (whether  or  not  it  was  con¬ 
firmed  by  the  Starovoit  et  al  bulletin  and  the  number  of  F.TDCs  that  reported  the  event).  We 
note  that  more  than  90%  of  the  entries  are  in  the  confirmed  category  (either  listed  by 
Starovoit  et  al  or  reported  by  at  least  one  EIDC). 

Table  7.63  summarizes  the  number  of  detected  events  by  the  various  systems.  We  note 
that  the  four  EIDCs  (reprocessed  bulletins  from  Stockholm,  Moscow,  Canberra  and 
Washington)  had  similar  performances,  and  reported  about  half  of  the  events  in  the  refer¬ 
ence  catalogue.  NE2C  reported  only  one  third  of  the  reference  events  in  their  monthly  bul¬ 
letin.  The  rapid  QED  service  (Quick  Epicenter  Determination)  reported  very  few  of  the 
events. 

The  GBF  association  process  reported  more  events  than  any  of  the  four  EIDCs,  and  also 
had  the  most  events  corresponding  to  the  reference  catalogue.  In  addition,  the  GBF 
method  produced  17  reports  that  did  not  correspond  to  entries  in  Starovoit  et  al’s  bulletin. 
Each  of  the  EIDCs  also  had  events  in  this  category,  but  not  as  many  as  the  GBF  process.  It 
should  be  noted  that  one  event  reported  by  one  EIDC  and  confirmed  by  Starovoit  et  al’s 
bulletin  was  not  reported  by  the  GBF  method.  The  reason  was  that  the  event  had  only  two 
valid  phases,  and  thus  did  not  satisfy  our  GBF  detection  criterion.  On  the  other  hand,  the 
GBF  reported  4  confirmed  events  that  were  not  in  any  of  the  EIDC  bulletins. 

We  also  conducted  an  experiment  to  test  the  likelihood  of  false  associations.  The  GBF 
process  with  the  parameters  used  in  this  study  was  run  on  a  7-day  period  prior  to  day  1 19. 
A  false  association  would  normally  correspond  to  phases  from  a  real  event  occurring 
somewhere  else,  but  for  which  the  phases  happened  to  match  our  criteria.  Table  7.6.4 
shows  the  events  associated  for  this  7-day  period.  Only  six  events  were  associated,  two  of 
which  were  in  fact  close  to  the  beam  steering  point.  Thus  only  four  definite  false  alarms 
were  observed  during  this  one-week  period.  We  conclude  that  the  false  alarm  rate  is  very 
low  for  this  processing  method. 

Conclusions 

The  GBF  technique  provides  a  simple  and  rapid  way  to  associate  large  numbers  of  phases 
from  an  aftershock  sequence  with  a  very  low  false  alarm  rate.  In  fact,  the  GBF  aftershock 
processing  of  24  hours  of  data  for  the  day  in  question  (29  April  1991)  took  only  5  minutes 
on  a  SUN  sparcstation2. 

We  consider  that  the  GBF  would  be  very  useful  as  a  preprocessor  to  the  expert  system 
algorithm  to  be  applied  at  a  future  International  Data  Center  (IDQ.  By  first  using  the  GBF 
to  extract  aftershock  sequences,  and  remove  the  corresponding  phase  detections,  the 
remaining  task  of  associating  events  from  other  locations  would  be  much  simplified.  Other 
applications  of  GBF  in  the  context  of  IDC  processing  can  also  be  envisaged.  Furthermore, 
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the  interaction  between  GBP  and  threshold  monitoring,  in  terms  of  eliminating  “unlikely” 
phase  associations,  deserves  to  be  studied  in  ''•tail. 

F.  Ringdai 
T.  Kvaerna 
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Station 

Type 

Lat 

Loo 

Distance 

oaa - 

nuw 

Trtime 

Azimuth 

Slowness 

KJV 

Single 

43.93 

42.68 

28.87 

157.33 

13.73 

KIV 

Single 

43.95 

42.68 

50.21 

157.33 

24.17 

ARU 

Single 

36.40 

58.60 

16.47 

p 

236.95 

221.31 

12.66 

GAR 

Single 

39.00 

70.30 

20.57 

p 

279.67 

288.38 

11.01 

FIN 

Hfarray 

61.44 

26.07 

21.66 

p 

291.44 

143.15 

10.52 

GER 

Hf array 

48.83 

13.70 

21.68 

tm 

291.58 

95.77 

10.51 

OSS 

Single 

46.69 

10.13 

24.02 

p 

315.08 

87.84 

9.61 

'ia 

Spanay 

60.13 

13.68 

25.33 

p 

327.43 

120.73 

9.28 

NRS 

Hf  array 

60.73 

11-54 

26.55 

p 

338.63 

118.86 

9.07 

ARC 

Hfarray 

69.54 

25.51 

28.66 

p 

353.11 

151.55 

8.92 

GBA 

Sparray 

13.62 

77.59 

41.16 

p 

465.27 

320.96 

8.23 

YKA 

Spanay 

62.49 

-114.61 

73.89 

p 

696.14 

16.68 

5.87 

Tkble  7.6.1.  Station  and  phase  parameters  used  for  GBF  processing  of  the  Caucasus  after¬ 
shock  sequence  (42.5N,  43.5E,  Depth  0). 
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Origin  tias 

Lit 

Lon 

Dspth 

Nph 

Mat* 

Tim 

Nssi 

Jurat 

Nslow 

Slraa 

Nslv 

Slvrss 

St arov 

EIDCS 

1991-119:09.  12.49.0 

42.50 

43.50 

0.00 

11 

11 

1.45 

7 

4.53 

7 

1.15 

7 

1.45 

Yss 

4 

1991-119:09.27.48.0 

42.50 

43.50 

0.00 

3 

3 

2.84 

2 

4.10 

2 

1.71 

2 

1.87 

Tat 

3 

1**1-119:09. 3X. 05.0 

42.50 

43.50 

0.00 

3 

3 

1.12 

3 

7.45 

3 

1.09 

3 

1.74 

NO 

3 

1*91-11*: 0*. 37 .39.0 

42.50 

43.50 

0.00 

11 

10 

2.29 

7 

4.88 

7 

1.24 

7 

1.64 

YSS 

2 

19*1-11*: 09. 38. 08.0 

42.50 

43.50 

0.00 

5 

5 

3.53 

4 

4.48 

3 

0.93 

3 

1.17 

No 

4 

1991-119:09.38.34.0 

42.50 

43.50 

0.00 

4 

4 

2.32 

4 

3.23 

3 

1.92 

3 

2.00 

NO 

3 

1991-119:09.41.52.0 

42.50 

43.50 

0.00 

4 

3 

1.23 

2 

7.25 

1 

0.93 

1 

2.02 

YSS 

3 

1*91-119:09.50.49.0 

42.50 

43.50 

0.00 

3 

2 

1.22 

1 

4.27 

0 

0.00 

0 

0.00 

NO 

2 

1991-119:09.54.37.0 

42.50 

43.50 

0.00 

8 

5 

2.87 

4 

5.33 

3 

1.01 

3 

1.42 

Yat 

4 

1991-119:09.59.24.0 

42.50 

43.50 

0.00 

11 

10 

1.12 

8 

4.43 

6 

0.90 

6 

1.19 

Yas 

4 

1991-119:10.01.15.0 

42.50 

43.50 

0.00 

9 

8 

2.32 

5 

2.94 

5 

0.94 

5 

1.05 

Yas 

4 

1991-119:10.08.23.0 

42.50 

43.50 

0.00 

4 

3 

3.93 

2 

2.84 

2 

1.70 

2 

1.77 

YSS 

3 

1991-119:10.08.37.0 

42.50 

43.50 

0.00 

3 

2 

1.02 

1 

4.27 

0 

0.00 

0 

0.00 

YSS 

_ 

1991-119:10.15.35.0 

42.50 

43.50 

0.00 

10 

9 

1.84 

5 

5.72 

4 

1.18 

4 

1.65 

YSS 

4 

1991-119:10.15.57.0 

42.50 

43.50 

0.00 

5 

5 

2.28 

4 

4.98 

3 

0.98 

3 

1.35 

No 

4 

1991-119:10.19.42.0 

42.50 

43.50 

0.00 

9 

8 

1.52 

5 

4.44 

5 

0.93 

5 

1.24 

YSS 

4 

1991-119:10.30.42.0 

42.50 

43.50 

0.00 

8 

5 

1.21 

4 

8.01 

3 

0.91 

3 

1.54 

YSS 

4 

1991-119:10.35  33.0 

42.50 

43.50 

0.00 

5 

4 

8.83 

2 

2.48 

1 

0.42 

1 

0.43 

YSS 

3 

1991-119:10.41.00.0 

42.50 

43.50 

0.00 

4 

3 

1.35 

2 

3.48 

1 

0.43 

1 

0.52 

YSS 

3 

1991-119:10.52.43.0 

42.50 

43.50 

0.00 

11 

10 

2.03 

8 

2.47 

6 

0.85 

6 

1.05 

YSS 

4 

1991-119:10.53.05.0 

42.50 

43.50 

0.00 

3 

3 

4.25 

2 

10.18 

2 

1.30 

2 

2.19 

NO 

_ 

1991-119:10.58.12.0 

42.50 

43.50 

0.00 

5 

4 

0.90 

3 

4.79 

2 

1.10 

2 

1.69 

YSS 

3 

1991-119:11.04.31.0 

42.50 

43.50 

0.00 

9 

8 

2.01 

4 

4.34 

4 

0.60 

4 

0.94 

YSS 

4 

1991-119:11.08.04.0 

42.50 

43.50 

0.00 

3 

3 

1.22 

3 

2.21 

2 

1.30 

2 

1.32 

NO 

2 

1991-119:11.10.14.0 

42.50 

43.50 

0.00 

10 

9 

2.34 

5 

8.35 

4 

0.72 

4 

1.42 

Yss 

4 

1991-119:11.12.21.0 

42.50 

43.50 

0.00 

3 

3 

0.39 

2 

7.52 

1 

0.94 

1 

2.26 

YSS 

3 

1991-119:11.39.38.0 

42.50 

43.50 

0.00 

4 

3 

0.48 

2 

2.29 

1 

0.14 

1 

0.14 

YSS 

3 

1991-119:11.43.19.0 

42.50 

43.50 

0.00 

5 

4 

2.02 

2 

2.48 

2 

0.66 

2 

0.81 

YSS 

4 

1991-119:11.5*  13 

42.50 

43.50 

0.00 

9 

8 

2.41 

5 

3.53 

4 

1.04 

4 

1.28 

Yss 

4 

1991-119:11. :>.r 6.0 

42.50 

43.50 

0.00 

10 

9 

2.98 

8 

8.76 

5 

0.77 

5 

1.57 

YSS 

4 

19*1-119:13.02 . 12.0 

42.50 

43.50 

0.00 

4 

3 

1.13 

2 

3.29 

2 

0.41 

2 

0.84 

YSS 

4 

1991-119:1'* .  12.39.0 

42.50 

43.50 

0.00 

6 

5 

1.79 

3 

5.06 

3 

0.23 

3 

0.97 

Yss 

4 

1991-119:^.13.26.0 

42.50 

43.50 

0.00 

4 

3 

2.05 

2 

6.46 

2 

1.01 

2 

1.50 

No 

4 

1991-119:13.19.50.0 

42.50 

43.50 

0.00 

8 

5 

0.80 

2 

2.29 

1 

0.14 

1 

0.11 

Yss 

4 

1?" 1-119: 13. 27. 17.0 

42.50 

43.50 

0.00 

9 

8 

2.84 

4 

7.83 

3 

0.33 

3 

1.67 

Yss 

4 

1991-119:13.49.59.0 

42.50 

43.50 

0.00 

7 

8 

0.87 

3 

2.94 

2 

0.66 

2 

0.73 

Yss 

4 

1991-119:13.53.10.0 

42.50 

43.50 

0.00 

5 

4 

0.82 

2 

2.93 

2 

0.72 

2 

1.03 

YSS 

4 

1991-119:14.00.28.0 

42.50 

43.50 

0.00 

4 

3 

2.38 

1 

0.68 

1 

0.20 

1 

0.21 

YSS 

4 

1991-119:14.20.57.0 

42.50 

43.50 

0.00 

3 

2 

0.49 

1 

4.27 

1 

0.82 

1 

1.10 

NO 

2 

1991-119:14.43.09.0 

42.50 

43.50 

0.00 

11 

10 

2.34 

6 

2.84 

6 

0.80 

6 

1.03 

YSS 

4 

1991-119:14.43.30.0 

42.50 

43.50 

0.00 

3 

3 

4.90 

2 

6.21 

2 

0.31 

2 

1.11 

NO 

- 

1991-119:15.28.48.0 

42.50 

43.50 

0.00 

6 

5 

1.29 

4 

7.45 

4 

0.95 

4 

1.75 

YSS 

4 

1991-119:15.38.58.0 

42.50 

43.50 

0.00 

4 

3 

2.38 

1 

11.77 

1 

1.90 

1 

2.73 

YSS 

2 

1991-119:18.03.09.0 

42.50 

43.50 

0.00 

8 

8 

5.82 

3 

2.55 

3 

1.03 

3 

1.14 

YSS 

3 

1991-119:18.12.49.0 

42.50 

43.50 

0.00 

5 

4 

1.55 

2 

5.70 

2 

1.06 

2 

1.47 

YSS 

3 

1991-119:18.22.27.0 

42.50 

43.50 

0.00 

5 

4 

1.45 

2 

2.29 

1 

0.64 

1 

0.84 

YSS 

4 

Nph  -  Number  of  associated  phases 

Nsta  -  Number  of  stations 

Tres  -  Mean  absolute  time  residual 

Nazi  -  Number  of  azimuth  observations 

Azres  -  Mean  absolute  azimuth  residual 

Nslow  -  Number  of  slowness  observations 

Sires  -  Mean  absolute  slowness  residual 

Nslv  -  Number  of  horizontal  slowness  vector  observations 

Slvres  -  Mean  absolute  horizontal  slowness  vector  residual 

Starov  -  Event  confirmed  by  Starovoit  et  al  catalogue  (Yes/No) 

EIDCs  -  Number  of  confirming  EIDCs 


Table  7.6.2.  List  of  event  parameters  for  the  events  detected  on  the  generalized  beam 
stewed  to  42.5°N,  43.5°E  for  day  119  (29  April)  1991.  See  text  for  details.  (Page  1 
of  2) 
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1991-119: 16.41.43.0 

42.50 

43.50 

0.00 

10 

9 

1991-119:1*. 49. 59.0 

42.50 

43.50 

0.00 

8 

7 

1991-119:1*. 58. 51.0 

42.50 

43.50 

0.00 

7 

6 

1991-119:17.10.29.0 

42.50 

43.50 

0.00 

6 

5 

1991-119:17.20.40.0 

42.50 

43.50 

0.00 

3 

3 

1991-119:17.21.27.0 

42.50 

43.50 

0.00 

7 

6 

1991-119:17.34.43.0 

42.50 

43.50 

0.00 

4 

3 

1991-119:17.55.01.0 

42.50 

43.50 

0.00 

9 

8 

1991-119:18.14.44.0 

42.50 

43.50 

0.00 

5 

4 

1991-119:18.17.22.0 

42.50 

43.50 

0.00 

4 

3 

1991-119:18.23.18.0 

42.50 

43.50 

0.00 

10 

10 

1991-119:18.30.43.0 

42.50 

43.50 

0.00 

9 

9 

1991-119:18.51.37.0 

42.50 

43.50 

0.00 

4 

3 

1991-119:19.07.05.0 

42.50 

43.50 

0.00 

12 

11 

1991-119:19. 16.04.0 

42.50 

43.50 

0.00 

6 

5 

1991-119:19.19.58.0 

42.50 

43.50 

0.00 

10 

9 

1991-119:19.26.52.0 

42.50 

43.50 

0.00 

3 

2 

1991-119:19.44.56.0 

42.50 

43.50 

0.00 

9 

8 

1991-119:19.52.52.0 

42.50 

43.50 

0.00 

4 

4 

1991-119:20.01.42.0 

42.50 

43.50 

0.00 

3 

2 

1991-119:20.12.08.0 

42.50 

43.50 

0.00 

9 

8 

1991-119:20.19.47.0 

42.50 

43.50 

0.00 

5 

4 

1991-119:20.24.45.0 

42.50 

43.50 

0.00 

10 

9 

1991-119:20.32.54.0 

42.50 

43.50 

0.00 

9 

9 

1991-119:21.23.16.0 

42.50 

43.50 

0.00 

4 

3 

1991-119:21.24.11.0 

42.50 

43.50 

0.00 

e 

8 

1991-119:21.25.24.0 

42.50 

43.50 

0.00 

4 

4 

1991-119:21.30.32.0 

42.50 

43.50 

0.00 

3 

2 

1991-119:22.25.07.0 

42.50 

43.50 

0.00 

3 

2 

1991-119:22.28.25.0 

42.50 

43.50 

0.00 

10 

9 

1991-119:23.10.54.0 

42.50 

43.50 

0.00 

4 

3 

1991-119:23. 17.56.0 

42.50 

43.50 

0.00 

3 

2 

1991-119:23.32.32.0 

42.50 

43.50 

0.00 

9 

8 

1991-119:23.32.50.0 

42.50 

43.50 

0.00 

3 

3 

1991-119:23.34.18.0 

42.50 

43.50 

0.00 

3 

3 

2.73 

5 

2.11 

5 

1.20 

5 

1.34 

1M 

4 

2.03 

3 

7.88 

2 

1.59 

2 

2.58 

TM 

4 

2.30 

3 

4.27 

2 

1.52 

2 

1.67 

YM 

4 

1.42 

2 

2.98 

1 

0.43 

1 

0.47 

Tee 

4 

0.57 

3 

6.21 

2 

0.89 

2 

1.55 

HO 

4 

1.99 

4 

3.32 

3 

1.48 

3 

1.63 

Tee 

4 

2.01 

1 

3.77 

1 

1.83 

1 

1.94 

Tee 

3 

2.07 

4 

3.45 

3 

1.17 

3 

1.28 

Tee 

4 

3.85 

2 

5.75 

1 

1.19 

1 

1.83 

Tee 

4 

2.92 

1 

1.73 

1 

1.05 

1 

1.09 

Tee 

2 

2.62 

6 

4.09 

6 

1.11 

6 

1.36 

Tee 

4 

2.37 

5 

3.86 

4 

1.22 

4 

1.38 

Tee 

4 

1.88 

2 

5.78 

2 

0.65 

2 

1.36 

Tee 

4 

1.59 

7 

3.81 

6 

1.17 

6 

1.38 

Tee 

4 

1.91 

2 

4.96 

2 

0.74 

2 

1.06 

Tee 

4 

2.78 

5 

4.10 

4 

0.77 

4 

1.22 

Tee 

4 

0.18 

0 

0.00 

0 

0.00 

0 

0.00 

Tee 

_ 

2.26 

4 

1.85 

4 

0.50 

4 

0.69 

Tee 

4 

1.40 

3 

2.24 

2 

0.18 

2 

0.24 

Tee 

4 

0.21 

1 

7.85 

1 

2.17 

1 

2.52 

Tee 

1.77 

4 

5.59 

4 

0.84 

4 

1.33 

Mo 

4 

1.91 

3 

3.06 

3 

0.99 

3 

1.16 

Tee 

4 

1.88 

5 

4.13 

4 

1.07 

4 

1.34 

Tee 

4 

0.87 

5 

3.50 

4 

1.48 

4 

1.63 

Tee 

4 

1.93 

1 

1.27 

1 

0.41 

1 

0.46 

MO 

2 

1.34 

5 

7.72 

4 

1.07 

4 

1.91 

Tee 

4 

0.95 

3 

8.45 

2 

0.77 

2 

1.96 

MO 

4 

2.58 

0 

0.00 

0 

0.00 

0 

0.00 

Mo 

3 

1.54 

0 

0.00 

0 

0.00 

0 

0.00 

Mo 

2 

2.62 

5 

5.00 

4 

0.61 

4 

1.29 

Tee 

4 

0.70 

1 

4.27 

0 

0.00 

0 

0.00 

Tee 

2 

1.13 

1 

4.27 

0 

0.00 

0 

0.00 

Tee 

- 

2.68 

4 

5.26 

3 

0.87 

3 

1.27 

Tee 

4 

6.58 

2 

4.89 

2 

0.70 

2 

1.27 

MO 

_ 

0.55 

2 

7.25 

1 

0.03 

1 

1.87 

Mo 

4 

Ibble  7.62  (amt).  (Page  2  of  2) 
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Source 

Total  number  of 
events 

Confirmed  by 
Starovoit  et  al’s 
catalogue 

Not  in  Starovoit  et 
al’s  catalogue 

Starovoit  et  al 
catalogue 

115 

115 

0 

Canberra  EIDC 
(reprocessed) 

57 

48 

9 

Stockholm  EIDC 
(reprocessed 

73 

62 

11 

Moscow  EIDC 
(reprocessed) 

76 

61 

15 

Washington  EIDC 
(reprocessed) 

71 

58 

13 

GBF  (automatic) 

82 

65 

17 

NEIC  monthly  list 

35 

35 

0 

QEDlist 

6 

6 

0 

Table  7.6.3.  Number  of  events  reported  by  various  sources  for  C s  W.  Caucasus  sequence 
of  29  April  1991.  From  our  analysis,  all  reported  GBF  events  for  that  day  were  real 
(no  false  alarms).  A  few  events  reported  by  the  BlDCs  or  GBF  were  close  in  time 
(possibly  multiple  events)  and  therefore  not  included  as  separate  events  in  Starovoit 
et  al’s  catalogue. 


a)  GBF  detection  list,  day  112-118 


NORSAR  Sci  Rep.  1-93/94 


the  six  events  were  actually  in  the  Caucasus  area,  while  the  four  remaining  detections  were  “side  lobes”  from  large  events  else¬ 
where  in  Eurasia.  Consequently,  there  were  only  4  “false  alarms”  for  the  entire  7-day  period. 
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Fig.  7.6.1,  Map  showing  the  stations  used  for  GBF  processing  of  the  Caucasus  aftershock 
sequence. 


